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INTRODUCTION 


The  qualitative  behavior  of  a  finite-length  periodic  array  of  small  antennas  or  scatterers 
closely  spaced  along  a  straight  line  can  be  understood  in  terms  of  traveling  waves,  supported 
by  the  corresponding  infinitely  long  periodic  array,  and  the  fields  of  the  feed  antenna  ele¬ 
ments)  of  the  array  that  are  not  converted  to  traveling  waves  [l]-[7].  Consider,  for  example, 
the  linear  array  (shown  in  Figure  1)  consisting  of  eight  perfectly  conducting  “dipoles”  each 
separated  by  a  distance  d  along  the  z  axis.  The  first  dipole  is  fed  by  a  sinusoidal  wave  with 
frequency  u>  and  the  other  seven  dipoles  axe  shorted  to  form  passive  straight-wire  scatterers. 
The  feed  dipole  radiates  time-harmonic  (e~wt,u  >  0)  fields  that  are  partially  converted  to 
a  traveling  wave  with  et/3z  dependence  for  discrete  values  of  z  along  the  array  separated  by 
integer  multiples  of  the  interelement  spacing  d.  (To  simplify  the  discussion,  assume  a  single 
traveling  wave  exists  at  the  frequency  u>.)  The  traveling  wave  propagates  without  decay  in 
the  positive  z  direction  with  the  real  propagation  constant  0  >  0  that  depends  on  frequency. 
If  the  electrical  size  of  each  element  is  small  enough,  the  periodic  array  of  elements  can 
support  only  a  slow  traveling  wave  with  phase  velocity  less  than  or  equal  to  the  speed  of 
light  such  that  0  >  k,  where  k  =  u>/c  =  2n/\  is  the  free-space  propagation  constant,  A  is  the 
free-space  wavelength,  and  c  is  the  speed  of  light.  Without  loss  of  generality,  we  can  assume 
that  0d  <  n  because  as  (3d  becomes  greater  than  n  the  traveling  wave  simply  becomes  a 
slow  wave  traveling  in  the  opposite  direction.1  Therefore,  (3d  satisfies  the  inequalities 

kd  <  (3d  <  7r  (1) 

which  imply  that  a  traveling  wave  can  exist  only  if  the  separation  distance  between  the 
electrically  small  elements  is  less  than  or  equal  to  half  the  traveling- wave  wavelength,  which 
is  always  less  than  or  equal  to  free-space  wavelength  A;  therefore  (d  <  A/2).  A  proof  of  the 
inequalities  in  (1)  for  a  general  linear  periodic  array  of  lossless  passive  small  elements  is  given 
in  Appendix  A  where  it  is  also  shown  that  the  fields  of  a  nonradiating  traveling  wave  on  a 
linear  infinite  periodic  array  must  decay  exponentially  with  radial  distance  from  the  axis  of 
the  array. 

The  traveling  wave  excited  by  the  feed  element  in  Figure  1  propagates  in  the  positive  z 
direction  and  diffracts  from  the  right  end  of  the  array  much  like  a  surface  wave  (or  waveguide 
mode)  diffracts  from  the  end  of  a  dielectric  rod  (or  open-ended  waveguide)  [8],  [9].  Part  of 
the  traveling  wave  is  reflected  and  part  is  radiated  unconfined  by  the  array.  The  reflected 
traveling  wave,  which  has  e~t)3z  dependence,  travels  the  length  of  the  array  and  diffracts  from 
the  left  end  of  the  array  in  Figure  1.  Again  part  of  the  traveling  wave  is  reflected  and  part 
is  radiated.  This  process  of  reflection  and  radiation  from  the  ends  of  the  array  continues  ad 
infinitum. 

If  the  propagation  constant  (3  of  the  traveling  wave  is  appreciably  greater  than  the  free- 
space  propagation  constant  k,  the  traveling  wave  is  strongly  attached  to  the  array  elements, 
the  impedance  mismatch  of  the  traveling  wave  as  it  encounters  free  space  at  the  ends  of  the 
array  will  be  large,  and  an  appreciable  amount  of  power  in  the  incident  traveling  wave  will 

xIf  the  elements  of  the  array  are  highly  dispersive,  it  is  possible  for  the  direction  of  power  flow  in  the 
traveling  wave  to  remain  in  the  +z  direction  for  0  >  nd.  If  we  then  change  this  positive  0  to  —0'  (for 
example,  if  ir/d  <  0  <  2n/d  let  0'  —  2n/d  —  0  <  ir/d),  the  traveling  wave  is  called  a  ’’backward  traveling 
wave”  [6,  pp.  264,  341-349], 
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Figure  1:  Linear  periodic  array  of  eight  perfectly  conducting  dipoles  separated  by  a  distance 
d  along  the  z  axis  and  fed  by  a  sinusoidal  source  with  frequency  u>. 
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be  reflected.  The  array  will  then  support  a  large  standing  wave  at  each  kd  value  for  which 
the  phase  of  the  traveling  wave  changes  by  an  integer  multiple  of  360°  as  the  traveling  wave 
starts  at  one  end  of  the  array,  travels  the  length  of  the  array,  reflects  from  the  other  end 
of  the  array,  returns  back,  and  reflects  from  the  starting  end  of  the  array.  The  traveling 
wave  will  interfere  constructively  near  these  kd  values  and  destructively  well  away  from 
these  values.  Therefore,  resonances  will  appear  in  the  plots  of  total  power  radiated  by  the 
array  (or,  equivalently,  total  power  accepted  by  the  feed  element  of  the  array)  versus  kd. 
As  the  value  of  /?  gets  closer  to  the  value  of  k,  the  traveling  wave  becomes  more  weakly 
attached  to  the  array  elements  and  is  launched  more  efficiently  from  the  ends  of  the  array  so 
that  the  resonances  caused  by  the  constructive  and  destructive  interference  of  the  traveling 
wave  become  less  pronounced.  Moreover,  as  the  value  of  /3  becomes  closer  to  that  of  k,  the 
significant  fields  of  the  traveling  wave  extend  an  increasingly  further  radial  distance  from 
the  axis  of  the  array.  Therefore,  the  effective  endfire  aperture  and  endfire  directivity  of  the 
array  will  increase  as  (3  gets  closer  in  value  to  k  (similar  to  the  increase  in  effective  aperture 
and  directivity  of  the  HEn  surface  wave  launched  from  the  end  of  a  dielectric  rod  [8]),  until 
/?  gets  so  close  in  value  to  k  and  the  traveling  wave  so  broad  that  it  cannot  be  efficiently 
excited  by  feeding  one  or  more  of  the  elements.  If  each  element  of  the  array  is  a  small 
resonant  antenna  (or  scatterer)  with  a  narrow  frequency  band  when  it  radiates  (or  scatters) 
alone  in  free  space,  this  single  element  resonance  will  be  superposed  on  the  resonances  of  the 
traveling  wave  [10]. 

The  resonances  of  the  traveling  wave  will  disappear  for  kd  greater  than  the  value  at  which 
the  amplitude  and  phase  of  the  coupling  between  the  array  elements  will  no  longer  support  a 
traveling  wave.  (According  to  (1),  this  value  is  always  less  than  or  equal  to  7r  for  small  array 
elements.)  For  kd  greater  than  this  value,  the  induced  excitation  of  each  passive  element 
diminishes  rapidly  with  increasing  distance  from  the  feed  element  to  the  passive  element, 
and  the  periodic  array  no  longer  radiates  as  an  “endfire”  antenna. 

The  foregoing  qualitative  description  of  the  operation  of  linear  periodic  arrays  can  be 
partly  gleaned  and  partly  surmised  from  the  experimental,  theoretical,  and  numerical  in¬ 
vestigations  into  dipole  (Yagi-Uda)  arrays,  as  exemplified  in  the  references  [1]— [7] .  It  is  not 
surprising  that  these  previous  investigations  were  limited  to  dipole  arrays  in  which  the  num¬ 
ber  of  dipoles,  their  spacing,  and  their  dimensions  did  not  vary  over  a  wide  range.  The 
accurate  construction  of  a  variety  of  linear  arrays  and  the  accurate  measurement  of  their 
near  and  far  fields  require  an  extremely  tedious  and  time  consuming  experimental  effort  [1]. 
Purely  theoretical  analyses  of  linear  arrays  have  been  limited  to  obtaining  approximate  so¬ 
lutions  to  integral  equations  for  infinite  dipole  arrays  [2],  [11].  Although  numerical  solutions 
to  approximate  integral  equations  have  been  successful  for  analyzing  finite  arrays  [3] ,  [4] , 
[5],  they  become  unwieldy  for  determining  the  effects  of  unrestricted  changes  in  the  array 
parameters,  especially  in  large  arrays.  Even  for  a  small  array,  an  accurate  numerical  solution 
to  the  integral  equations  is  very  difficult  to  obtain  if  the  individual  elements  of  the  array 
are  resonant  antennas  or  scatterers  [10].  Moreover,  the  numerical  difficulties  increase  if  the 
elements  of  the  array  do  not  lie  in  a  straight  line,  such  as  in  closed-loop  arrays  [12],  or  if  the 
elements  are  distributed  over  a  surface  or  throughout  a  volume. 

The  scattering-matrix  analysis  of  linear  periodic  arrays  of  electric  dipoles  given  in  this 
report  was  suggested  by  the  corresponding  analysis  performed  by  Yaghjian  [13]  for  linear 
periodic  arrays  of  electroacoustic  monopole  transducers.  The  basic  ideas  of  this  report  were 
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worked  out  by  Yaghjian  within  the  simpler  analytic  framework  of  scalar  spherical  wave 
functions.  The  extension  to  electromagnetic  antennas  reported  here  requires  the  use  of  the 
more  complicated  vector  spherical  wave  functions.  However,  no  new  concepts  are  needed 
for  our  treatment  of  electromagnetic  antennas.  In  this  report  we  formulate  a  general  vector 
spherical-wave  scattering-matrix  representation  of  antennas  and  derive  reciprocity  and  loss¬ 
less  power  conservation  relations  between  the  antenna  reflection,  receiving,  transmitting,  and 
scattering  coefficients  that  constitute  the  scattering  matrix.  We  then  assume  that  each  array 
element  is  an  electrically  small  dipole  antenna.  This  specialization  of  the  general  framework 
to  small  electric  dipole  antennas  considerably  simplifies  the  source  scattering-matrix  analysis 
since  only  two  vector  spherical  wave  modes  (one  incoming  and  one  outgoing)  are  required  to 
describe  the  antennas.  The  scattering  matrix  involves  the  receiving,  transmitting,  and  scat¬ 
tering  coefficients  of  each  antenna  element,  as  well  as  the  input  reflection  coefficient  of  the 
waveguide  mode  feeding  each  antenna.  However,  for  a  linear  periodic  array,  the  reciprocity 
and  lossless  power  conservation  relations  reduce  the  critical  scattering-matrix  variables  in 
the  TV  equations  for  the  TV  unknown  outgoing  wave  coefficients  of  the  array  elements  to  one 
single  parameter,  the  phase  of  the  effective  scattering  coefficient  of  each  antenna  element. 
(The  number  of  elements  TV  and  the  normalized  spacing  kd  of  the  elements  in  the  linear 
periodic  array  are  the  other  critical  variables  in  the  equations.) 

The  magnitude  of  this  effective  scattering  coefficient  does  not  enter  the  equations  because 
reciprocity  and  power  conservation  demand  that  the  magnitude  |5e|  and  phase  V’e  of  the 
effective  scattering  coefficient  (Se  =  |5e|e‘^e)  satisfy  the  simple  relationship 

|5e|  =  ^sin^e  (2) 

which  implies  that  0  <  ipe  <  ir  and  that  ipe  approaches  0  or  7r  as  the  magnitude  of  the 
scattering  coefficient  approaches  zero.  We  also  prove  regardless  of  reciprocity  and  power 
conservation  relations  that  an  infinite  linear  periodic  array  of  short  electric  dipole  antennas 
oriented  perpendicular  to  the  array  axis  can  support  a  nondecaying  traveling  wave  only 
if  the  relationship  (2)  is  satisfied.  Remarkably,  in  contrast  with  the  acoustic  monopole 
case,  it  is  possible  to  have  two  traveling  waves  for  each  value  of  ^e.  A  simple  implicit 
transcendental  equation  determines  the  kd  —  (3d  diagram  for  these  traveling  waves.  The 
theoretical  propagation  velocity  of  the  traveling  wave  supported  by  a  periodic  array  of  short 
thin  wires  obtained  from  the  kd  —  (3d  diagram  in  combination  with  values  of  the  effective 
scattering  coefficient  of  the  wires  calculated  using  the  NEC  computer  code  [26]  are  shown  to 
compare  well  with  the  measurements  of  Ehrenspeck  and  Poehler  [1]  for  the  case  of  Yagi-Uda 
arrays. 

All  of  the  properties  exhibited  by  finite  linear  periodic  arrays  explained  in  the  first  part 
of  this  introduction  are  clearly  and  conveniently  demonstrated  by  solving  the  TV  x  TV  matrix 
equation  for  these  arrays  as  the  value  of  the  phase  V’e  of  the  effective  scattering  coefficient  is 
varied  from  0  to  it  and  the  value  of  kd  ranges  from  0  to  n  (that  is,  over  the  range  of  kd  at 
which  a  traveling  wave  can  exist).  In  addition,  we  confirm  the  result,  found  experimentally 
by  Ehrenspeck  and  Poehler  [1]  and  predicted  theoretically  by  Hansen  and  Woodyard  [14], 
that  the  maximum  endfire  directivity  attainable  by  varying  the  number  of  elements  in  a 
linear  periodic  array  of  electrically  small  elements  supporting  a  single  traveling  wave  (and 
having  one  feed  element)  depends  only  on  (and  increases  monotonically  with)  the  relative 


4 


phase  velocity  of  the  traveling  wave. 

The  remainder  of  the  paper  is  divided  into  eight  main  sections  in  addition  to  the  Conclud¬ 
ing  Remarks  (Section  10)  and  the  Appendices.  The  reader  interested  in  just  the  equations 
and  results  for  linear  periodic  arrays  of  small  electric  dipoles  oriented  perpendicular  to  the 
array  axis,  and  not  in  how  they  are  derived,  can  concentrate  on  Sections  6  and  7,  which  dis¬ 
cuss  infinite  and  finite  arrays,  respectively.  Section  2  introduces  the  basic  field  equations  and 
the  vector  spherical-wave  source  scattering-matrix  description  of  a  general  electromagnetic 
antenna  system.  Sections  3  and  4  derive  reciprocity  and  lossless  power  relations,  respec¬ 
tively,  that  must  be  satisfied  by  the  coefficients  of  the  scattering  matrix  of  a  general  linear, 
reciprocal,  lossless  antenna.  Section  5  applies  these  reciprocity  and  power  relations  to  the 
small  electric  dipoles  that  are  used  in  Sections  6  and  7  for  the  analysis  of  infinite  and  finite 
arrays  of  dipoles  oriented  perpendicular  to  the  array  axis.  Sections  8  and  9  consider  arrays 
of  dipoles  aligned  with,  and  skew  to,  the  array  axis. 


2  SOURCE  SCATTERING-MATRIX  DESCRIPTION 
OF  A  GENERAL  ANTENNA 

The  electromagnetic  system  under  consideration  is  pictured  in  Figure  2.  The  antenna  is 
bounded  by  the  closed  surface  Ap.  The  electromagnetic  power  supply  or  detector  is  bounded 
by  the  closed  surface  Ap  and  the  power  supply  feeds  the  antenna  through  a  waveguide  with 
reference  plane  A0  common  to  both  At  and  Ap.  The  power  supply  sends  time-harmonic 
(e_ta,<,w  >  0)  electromagnetic  power  down  the  waveguide  in  the  form  of  a  single  propagating 
waveguide  mode.  (The  waveguide  is  assumed  to  be  composed  of  perfect  conductors  separated 
by  a  linear,  homogeneous,  isotropic  medium.)  The  antenna  converts  the  electromagnetic 
fields  of  the  waveguide  mode  into  time-harmonic  electromagnetic  fields  that  propagate  in 
the  free  space  in  which  the  antenna  and  its  power  supply  or  detector  are  placed.  The  space 
between  the  spherical  reference  surface  Ar  with  radius  r  and  the  external  surface  Ap  =  (Ap— 
A0)  +  (Ap  —  A0)  of  the  antenna  plus  power  supply  or  detector  contains  no  electromagnetic 
sources,  although  arbitrary  time-harmonic  electromagnetic  sources  at  frequency  u>,  such  as 
other  antennas,  may  exist  outside  Ar.  It  is  assumed  that  the  external  part  of  the  surface  of  the 
power  supply  or  detector  (Ap—A0)  is  electromagnetically  shielded  so  that  no  electromagnetic 
power  crosses  (Ap  —  Aq). 

The  Maxwell  equations  that  govern  the  time-harmonic  electromagnetic  fields  in  the 
source-free  free  space  between  the  surface  Ap  and  the  spherical  surface  Ar  are  [15,  sec. 
7.1] 

V  x  E(r)  =  ikZoH  (3a) 

V  x  H(r)  =  -ikE(r)/Z0  (3b) 

which  imply  that  for  k  ^  0 

V  •  E(r)  =  0  (3c) 

V  •  H(r)  =  0  (3d) 

where  E(r)  and  H(r)  are  the  electric  and  magnetic  fields,  Z0  is  the  free  space  impedance, 
k  =  ujjc  =  2tt/A  is  the  free  space  propagation  constant,  A  is  the  free-space  wavelength,  and 
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Figure  2:  Schematic  of  the  electromagnetic  antenna  system.  The  external  surface  of  the 
power  supply  or  detector  is  given  by  Ae  =  {At  —  A$)  +  (Ap  —  Ao). 


c  is  the  speed  of  light.  Eliminating  either  H(r)  or  E(r)  from  the  Maxwell  equations  (3a)  and 
(3c)  yields  the  vector  wave  equations 

-V  x  V  x  E(r)  +  k2  E(r)  =  0  (4a) 

— V  x  V  x  H(r)  +  k2  H(r)  =  0  (4b) 

which  imply 

V2E(r)  +  k2  E(r)  =  0,  V-E(r)=0  (5a) 

V2H(r)  +  k2  H(r)  =  0,  V-H(r)  =  0.  (5b) 

Since  the  electric  field  satisfies  (5a)  across  the  surface  Ar  (which  lies  within  free  space), 
the  electric  field  on  the  surface  can  be  expanded  in  a  complete  set  of  divergenceless  vector 
spherical  wave  functions  [15], [16] 


The  spherical-wave  expansion  of  the  magnetic  field  then  follows  from  (3a) 

OO  l 

H(r)  =  —  ir„E  E  [«2Nffi(r)  +  CN2(r)+«SMffi(r)+tSMg(r)  ,  (7) 

|=1  m=—l 

where  Vo  =  1/Z0  is  the  free-space  admittance.  The  vector  spherical  harmonics  M^(r)  and 
N^(r)  are  defined  in  Appendix  B.  The  M^(r)  spherical  wave  functions  have  no  radial 
component.  Because  of  that  the  M^(r)  functions  in  (6)  and  the  N^(r)  functions  in  (7)  are 
often  referred  to  as  TE  waves  while  the  N^(r)  functions  in  (6)  and  the  M^(r)  functions  in 
(7)  are  often  referred  to  as  TM  waves;  they  satisfy  the  equations 

N£(r)  =  jV  x  Mj«(r)  (8a) 

M£(r)  =  jV  x  Ng>(r).  (8b) 

The  (r)  and  (r)  have  the  radial  dependence  ji(kr)  while  the  (r)  and  N[^| (r)  have 
the  radial  dependence  h\x\kr)  where  j[(kr)  and  hj^(fcr)  are  the  spherical  Bessel  functions 
and  the  spherical  Hankel  functions  of  the  first  kind,  respectively.  The  expansions  in  (6)  and 
(7)  are  known  as  multipole  expansions.  The  l  —  l  terms  are  the  dipole  terms,  the  l  =  2  terms 
are  the  quadrupole  terms,  etc.  The  summations  containing  the  modal  coefficients  in  (6) 
and  (7)  equal  the  electromagnetic  fields  produced  by  the  sources  (applied  and  induced)  that 
reside  outside  the  surface  Ar.  The  summations  containing  the  modal  coefficients  in  (6) 
and  (7)  equal  the  electromagnetic  fields  produced  by  the  sources  (applied  and  induced)  that 
reside  inside  the  surface  Ar  or,  equivalently,  inside  the  external  surface  Ae  of  the  antenna 
plus  power  supply  or  detector.  They  are  the  fields  radiated  and  scattered  by  the  antenna 
and  its  power  supply  or  detector. 
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The  tangential  electric  and  magnetic  fields  (Etan(p),  Htan(p)]  on  the  reference  plane  Ao 
of  the  feed  waveguide  can  be  written  in  terms  of  the  electric  and  magnetic  basis  fields 
[e0(p),  ho(p)]  of  the  single  propagating  waveguide  mode  with  incident  amplitude  a0  and 
emergent  amplitude  60;  specifically 

E  tan(p)  —  K>e0(p),  Htan(p)  =  Zoho  (p)  (®) 

where 

Vo  =  o0  +  bo,  Io  =  %(°o  —  bo).  (10) 

There  may  be  evanescent  modes  in  the  feed  waveguide,  but  the  fields  of  these  evanescent 
modes  are  assumed  to  be  negligible  on  the  reference  plane  Aq. 

If  the  dimensional  units  of  e0  and  ho  are  chosen  as  (meter)-1  and  they  are  consistent 
with  Maxwell’s  equations  in  the  International  System  of  mksA  units,  then  Vo  has  units  of 
Volts,  Io  has  units  of  Amperes,  and  the  characteristic  admittance  r)0  of  the  waveguide  can  be 
chosen  as  a  real  positive  constant  with  units  of  (Ohm)-1.  It  then  follows  that  normalization 
of  the  basis  fields  may  be  expressed  as  a  nondimensional  number  equal  to  one,  that  is 

J  Mp)  x  Vp)]  ’  *iodA  =  1  (n) 

Ao 

where  is  the  inward  normal  to  the  antenna  on  the  plane  Ao-  If  the  plane  Ao  simply  cuts 
two  wire  leads  at  quasi-static  frequencies,  Vo  and  Io  refer  to  conventional  circuit  voltages 
and  currents  that  do  not  serve  as  genuine  modal  coefficients.  In  that  case,  the  equations  in 
(10)  become  definitions  of  a0  and  b0  such  that 

ao  =  (Vo  +  W%)/2>  b0  =  (Vo  —  Io/rj0)/2.  (12) 

The  “outgoing”  modal  coefficients  (bo,b\H)  of  the  antenna  system  are  related  to  the  “in¬ 
coming”  modal  coefficients  (a0,  a££)  of  the  antenna  system  by  a  linear  matrix  transformation 
termed  the  “source  scattering-matrix  equations”  [17],  [18]  for  the  spherical- wave  representa¬ 
tion  of  the  antenna.  These  source  scattering-matrix  equations  can  be  written  as 


b0=Ta0+Y:  £  £  *; 


l=\  m=—l  s=l 


(13) 


=  r,(s) 


lma 0 


i'= 1  m'=-V  s'= 1 


(14) 


The  coefficients  R^,  T^,  and  Sj£:jfr 2>  in  the  source  scattering  matrix  embody  the  receiving, 
transmitting,  and  scattering  properties  of  the  linear  antenna.  The  complex  number  T  is  the 
reflection  coefficient  (phase  referenced  to  the  plane  Ao)  of  the  waveguide  mode  and  equals 
b0/a0  when  all  sources  outside  the  external  surface  Ap  of  the  antenna  and  its  power  supply 
are  zero,  in  other  words,  when  there  are  no  sources  outside  Ar  so  that  all  the  afm  equal  zero. 

The  spherical-wave  source  scattering  matrix  given  in  (13)-(14)  for  electromagnetic  anten¬ 
nas  differs  from  the  classical  spherical-wave  scattering-matrix  given  in  [19,  secs.  9.18-9.24] 
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because  the  classical  scattering  matrix  uses  (h\2\h\^)  radial  basis  functions  in  (6)  and  (7) 
instead  of  (j(,  h\^).  This  small  mathematical  change  from  tif*  to  ji  appreciably  simplifies 
the  analysis  by  relating  all  the  scattering-matrix  coefficients  directly  to  the  physical  sources. 
For  example,  if  an  antenna  with  its  power  supply  or  detector  does  not  scatter  fields  produced 
by  external  sources,  then  the  source  scattering  coefficients  Sim-i'm'  8X8  equal  to  zero.  This 
is  not  the  case  for  the  classical  spherical-wave  scattering-matrix  description  of  antennas  [20, 
sec.  2.3.5]. 

Finally,  note  that  for  a  scatterer  rather  than  an  antenna,  we  can  set  and  equal 
to  zero  and  equation  (13)  becomes  superfluous. 


3  RECIPROCITY  RELATIONS 

The  receiving,  transmitting,  and  scattering  coefficients  in  the  source  scattering  matrix  for 
antennas  composed  entirely  of  reciprocal  material  must  satisfy  reciprocity  relations  that  can 
be  derived  from  the  following  electromagnetic  reciprocity  theorem  [21,  sect.  1.9] 

/  [E 2(p)  x  Hx(p)  -  Ei(p)  x  H 2(p)]  •  fio<M  =  /  [E2(r)  x  Hx(r)  -  Ex(r)  x  H2(r)]  •  MA 

Ao  Ar 

(15) 

where  [Ei(p),Hi(p),Ei(r),Hi(r)]  and  [E2(p),H2(p),E2(r),H2(r)]  are  the  fields  of  the  an¬ 
tenna  excited  by  any  two  sets  of  electromagnetic  sources  located  inside  the  power  supply  or 
outside  the  spherical  surface  Ar. 

Substituting  the  electric  and  magnetic  fields  on  A0  in  the  waveguide  from  (9)-(10)  into 
the  left-hand  side  (LHS)  of  (15),  and  then  using  the  integral  relation  (11)  we  find  that  the 
LHS  of  (15)  is  equal  to 

2%  (&20a i0  -  Ka2o)  •  (16) 

Substituting  the  electric  and  magnetic  fields  on  Ar  from  (6)  and  (7)  and  using  the  orthogo¬ 
nality  relations  [20,  Eq.  Al.71]  (see  (124), (125)) 


2tt  7T 

/  /  [m2M  x  Mg(r 
0  0 

)|  •  r  sin  QdOdcj)  =  0 

(17a) 

27T  7T 

/  /  [NSW  x  N»(r) 
0  0 

ij  •  r  sin  6d0d(f)  =  0 

(17b) 

we  see  that  all  integrals  of  cross-products  involving  Mj^(r)  x  M^(r)  and  N^(r)  x  N^](r) 
can  be  set  equal  to  zero.  If,  in  addition,  we  make  use  of  the  reciprocity  integral  relations  [20, 
Eqs.  A1.72-A1.74]  (see  (124), (125)) 

/  [Mg(r)  x  Ng(r)  -  Mg>(r)  x  N®(r)]  ■  <IS  =  (18a) 
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/  i\m+l 

/  [N<«(r)  x  M«(r)  -  N“(r)  x  M<«(r)]  •  dS  = 


(18b) 


where  A^1’1)  =  0,  A^1’2)  =  i,  A^2’1)  =  — i,  A^2,2^  =  0  and  the  surface  of  integration  in  (18a)  and 
(18b)  is  a  spherical  surface  S  of  radius  r  with  outward  area  element  dS  =  rr2  sin  dd9d<j>,  we 
find  that  the  right-hand  side  (RHS)  of  (15)  is  equal  to 

*0  V2'  V'  (  1  lm+1  ro(1)  hW  —  +  a ft®  —  bP  aP}  1  ( 191 

TJ  2-x  Z-x  V—l/  °2; '  a2;lm°l;l,-m  0'2;l,ma2;l,-rn\  Viy/ 

K  1=1  m=-J 

oo  l  2 

Letting  E  =  E  E  E  and  equating  (16)  with  (19)  we  obtain 

Ims  J=1  m=—l  5=1 

2%  (&2o«io  -  =  g  E(-i)"+I  [<45>&U  -  4'WtJ  ■  (M) 

Ims 

Substituting  expressions  for  620  >  &i0>  and  62*jm  fr°m  the  source  scattering-matrix  equa¬ 

tions  (13)  and  (14),  the  equation  (20)  becomes 


OloE  \2W 


*)  (-l)m+1^o 


t>(«)  1  _  flo  v  [277  r(s)  -  (  i)w+1y0/3)  1  ( ) 

*l-m  a2 ;lm  a2o  2-^  Z7iOnlm  £2  1l-m  al ;lm 

Ims  L 


V  „(s)  (  i)m+loW,(j')  -(s') 


ims  /'mV 


+§E  E  4tn(-ir+1^a<45L =o. 


Ims  /'mV 


By  rewriting  the  indices  of  the  fourth  term  we  then  obtain 


o„  r>(»)  (”1)  Y0m(s)  _(s)  _  d(«)  (  jO  „W 

Olo  X!/  ^Vo^lm  U2  °2;/m  °2o  ^  ^Vo^lm  jU2  -^l,-m  al ;/m 

Ims  L  *  J  ims  L  *  J 

-§  E  E  <41  [(-ir«s«S>„,  -  (-ir'+1s£2&J  =  o.  (22) 

^  ims  /'mV 


Since  the  modal  coefficients  ai0  and  a2o  are  determined  by  the  power  supply  connected  to  the 
waveguide,  and  the  modal  coefficients  a^m  and  a2®^  are  determined  by  the  sources  outside 
Ar,  they  can  all  be  chosen  independently.  Therefore  (22)  can  hold  for  all  values  of  the 
modal  coefficients  if  and  only  if  the  quantities  in  each  of  the  sets  of  square  brackets  are  zero. 
Specifically,  the  receiving,  transmitting,  and  scattering  coefficients  of  the  source  scattering 
matrix  must  satisfy  the  reciprocity  relations 

2-foRS  =  {-=% (23) 

^i,-m;/'m'  *”  \  l)  ^/',-m';Zm- 

The  reciprocity  relation  (23)  states  that  the  spherical-wave  receiving  coefficients  of  an 
antenna  containing  only  reciprocal  material  are  equal  to  a  constant  times  the  spherical- wave 


10 


transmitting  coefficients.  The  reciprocity  relation  (24)  states  that  the  fields  scattered  into 
one  outgoing  spherical  mode,  when  another  spherical  mode  is  incident  upon  the  antenna, 
remain  the  same  if  the  mode  numbers  of  the  incident  and  outgoing  modes  are  interchanged. 
For  a  scatterer  rather  than  an  antenna,  =  0  and  only  the  reciprocity  relation 

(24)  is  relevant. 


4  POWER  CONSERVATION  RELATIONS 

Conservation  of  power  also  places  restrictions  upon  the  coefficients  of  the  scattering  matrix. 
If  the  antennas  are  lossless,  the  net  average  power  flowing  across  the  reference  plane  Aq  of  the 
waveguide  toward  the  antenna  must  equal  the  net  average  power  flowing  out  of  the  spherical 
surface  Ar.  Mathematically,  this  conservation  of  average  power  flow  can  be  expressed  from 
Poynting’s  theorem  as 

^Re  J  (Etan(p)  x  HIM]  *  %£A  =  ^Re  J  [E(r)  x  H*(r)]  •  dS.  (25) 

Aq  Ar 

Inserting  Eton(p)  and  Htan(p)  from  (9)-(10)  into  the  left  hand  side  of  (25),  and  then  using 
the  integral  relation  (11)  the  LHS  of  (25)  is  found  to  equal 

(kf  -  M2)  (26) 

To  evaluate  the  RHS  of  (25)  we  begin  by  substituting  the  vector  spherical  wave  function 
expansions  of  E(r)  and  H(r)  given  by  (6)  and  (7).  We  then  use  the  orthogonality  relations 
[16,  Eqs.VI-7,VI-8]  together  with  [15,  Eq.9.121],  and  [16,  Eqs.VI-6,VI-7]  together  with  [15, 
Eq.9-120],  along  with  the  Wronskian  formulas  for  spherical  Bessel  and  Hankel  functions  to 
show  that  the  RHS  of  (25)  is  equal  to 

E  [1^  I2  +  R*  (<&£')]  •  (27) 

Equating  (26)  and  (27)  we  obtain  the  conservation  of  power  relation 

(kl2  -  IM2)  =  515^0  E  [l<fi!l2  +  Re  (<&£*)]  •  (28) 

^  *  Ims 

If  the  antenna  neither  radiates  nor  scatters,  the  bj^  modal  coefficients  would  all  be  zero  and 
the  net  average  power  crossing  the  sphere  Ar  would  be  zero,  even  though  the  modal 
coefficients  generated  by  electromagnetic  sources  outside  Ar  may  not  be  zero.  In  that  case, 
the  average  power  entering  the  surface  Ar  would  equal  the  average  power  leaving.  With  the 
right  hand  side  of  (28)  zero,  the  left  hand  side  must  also  be  zero  for  this  lossless  antenna, 
and  all  the  power  in  the  waveguide  mode  propagating  toward  the  antenna  is  reflected  back 
down  the  waveguide  (|&o|2  =  |o0|2)- 
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When  the  bo  and  b^l  modal  coefficients  from  the  source  scattering-matrix  equations  (13) 
and  (14)  are  substituted  into  (28),  this  equation  can  be  written  as 

\Vok2 


kl 


>o 


(i-|r|2)-E|itT 


Ims 


Vpk2 

Y0 


Ims 


+  E 

Ims  \\ 


l\  V-  q(sW)  (S') 

I  ^Im^Um^Um* 


+  Re  E  Ste 

Urn's* 


>.(*')„(«')  „(*)*] 

•J'm'Q'l'm'Q'lm  J 

2ReLE(^r^*oW*  +  ir«aS-+  E  itimfj’altt)]  =  0. 

.  Ims  \  r°  Z  I'm' s'  /J 


(29) 


Because  the  modal  coefficients  a0  and  a\^  can  be  chosen  independently,  let  all  the  equal 
zero  to  get  from  (29)  the  first  power  relation 

^  (1  -  in2) . 


E|3 


M 


Yev  -./■  f30) 

Ims  10 

Next,  let  a0  =  0  to  show  that  the  quantity  inside  the  second  set  of  square  brackets  in  (29)  is 
zero.  Since  all  but  any  one  of  the  aj^  can  be  chosen  zero,  the  quantity  in  the  second  set  of 
square  brackets  in  (29)  can  vanish  for  all  only  if 


Dp  f  q(*).(«)\  .  V  I  I2  —  |n(.)| 

Ke  )  +  2L,  ~  y  \nlm\ 


(31) 


I'm's' 


Lastly,  the  quantity  inside  the  third  set  of  square  brackets  in  (29)  must  be  zero  because  the 
quantities  inside  the  other  two  sets  of  square  brackets  are  zero.  All  but  one  a0a^2  a  ^me 
can  be  chosen  zero  and  the  one  nonzero  a0°!m  can  be  given  a  purely  real  or  purely  imaginary 
value  in  order  to  reveal  a  third  power  relation 


rps 
^  lm 


I  rp(s')  q(s'),(s)*  _ 

+  2-j  1l'm'°lm-,Vm'  ~ 
Urn's* 


•lo^PpW* 

_i  Klm 


Yo 


(32) 


The  power  relation  (30)  expresses  mathematically  that,  in  the  absence  of  exterior  sources, 
the  total  power  radiated  by  a  lossless  antenna  equals  the  power  supplied  to  the  antenna 
through  the  waveguide.  The  power  relation  (31)  implies  that  the  total  power  received  by  a 
lossless  antenna  that  is  not  transmitting  can  always  be  related  to  the  scattering  coefficients 
of  the  source  scattering  matrix.  The  power  relation  (32)  shows  that  for  lossless  antennas 
the  product  of  the  reflection  coefficient  and  each  receiving  coefficient  of  the  source  scattering 
matrix  is  linearly  related  to  the  corresponding  transmitting  coefficient  and  the  sum  of  prod¬ 
ucts  of  all  the  transmitting  and  scattering  coefficients.  This  last  power  relation,  which  comes 
from  the  cross  coupling  terms  in  the  quadratic  power  expression  (29),  has  no  immediately 
obvious  physical  interpretation. 

For  lossless  scatterers  rather  than  lossless  antennas,  =  0  and  |T|  =  1  so  that 

the  power  relations  (30)  and  (32)  vanish  and  (31)  becomes 

(33) 


Re  (sag)  +  E  =  0- 


Urn's* 


12 


SHORT  Z  DIRECTED  ELECTRIC  DIPOLES 


We  now  specialize  the  general  vector  spherical  wave  expansions  (6)  and  (7)  and  the  source 
scattering-matrix  equations  (13)  and  (14)  to  the  case  of  a  short  z  directed  electric  dipole 
antenna.  The  electric  and  magnetic  fields  of  an  elementary  electric  (Hertzian)  dipole  antenna 
are  proportional  to  (r)  and  — iloM^  (r),  respectively.  Explicit  expressions  for  these 
outgoing  vector  spherical  waves  are  [16,  pp.  40-41] 


(r)  = 


ih^(kr) 


cos  Or 


kr  d(kr) 


[Ax/i^/cr)]  sin 06, 


— iloMio*  (r)  =  Yoh?(kr)y~sm04>. 


(34b) 


Substituting  expressions  for  the  spherical  Hankel  functions,  we  find  the  r  and  0  components 
of  Njo^  (r)  are  given  by 

(35a) 

Nl»  - -£  (i  • +  i -  w?)  Its**  ‘35b> 


and  similarly 


-iYoMZ(r)  =  -Yo^r  (l  +  smO. 


Equations  (35a)- (35b)  and  (35c)  are  the  components  of  the  electric  and  magnetic  fields, 
respectively,  of  a  z  directed  electric  point-current  source  (Hertzian  dipole)  [22,  p.  436].  The 
electric  fields  incident  upon  the  antenna  from  sources  outside  Ar  are  given  by  the  incoming 
(r)  and  Mj^(r)  terms  in  (6).  If  the  dimensions  of  the  antenna  (and  its  power  supply  or 
detector)  are  appreciably  smaller  than  a  wavelength,  and  the  origin  of  the  coordinate  system 
is  centered  on  the  antenna,  the  incident  electric  field  over  the  region  of  space  occupied  by 
the  antenna  (r  <C  A)  can  be  well  approximated  by  just  the  NjJ^(r)  modes  because  the 
electric  fields  of  all  the  other  modes  approach  zero  as  kr  approaches  zero.  Additionally,  as 
kr  approaches  zero  N^±1(r)  have  only  an  x  and  y  electric-field  component.  Moreover,  if  we 
assume  that  the  scattered  fields  of  the  electric  dipole  antenna  are  also  electric  dipole  fields, 
(6)  shows  that  the  z  component  of  the  electric  field  just  outside  the  antenna  is  given  by 

Ez( r)  =  [aio  Nffi(r)  +  A^N^r)]  •  z,  kr  <  1  (36) 


which  can  be  expanded  as 


E,(r)  =  oS’ii 


27T  kr 


^  1  /  i 

“‘nrs 


fcr<  1.  (37) 
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If  we  let  E °  refer  to  the  incident  electric  field  at  the  center  (origin)  of  the  antenna,  that  is 


(38) 


and  let  6  refer  to  the  coefficient  of  the  elkr/ ( kr )  radiated  electric  field,  that  is 


(39) 


(37)  becomes 


=  I?  [*  {Ir  +  (W)  +  (l +  1  ~ 


(fcr)2 


sin  00 


>z,  kr  «  1 .  (40) 


We  note  that  the  6-term  in  (40)  that  is  dotted  into  z,  which  gives  the  electric  field  radiated 
and  scattered  by  the  small  electric  dipole  antenna,  holds  for  all  kr  outside  the  antenna. 

Since  all  the  are  zero  except  for  an  antenna  that  radiates  as  an  electric  dipole, 
the  reciprocity  relation  (23)  demands  that  all  the  are  zero  except  R^.  Also,  with  all 
the  6^2  and  zero  except  when  Z  =  l,m  =  0,  and  s  =  2,  (14)  and  the  reciprocity  relation 
(24)  imply  that  all  the  Si'm'M  are  zero  except  S$i(02)-  Thus,  the  source  scattering-matrix 
equations  (13)  and  (14)  for  such  an  electrically  small  z  directed  dipole  reduce  to 


b0  =  Ta0  +  RE°Z  (41) 

b  =  Ta0  +  SE°z  (42) 

in  which  the  renormalized  receiving,  transmitting,  and  scattering  coefficients  for  the  small  z 
directed  electric  dipole  antenna  are  defined  as 


The  reciprocity  relation  (24)  reduces  to  an  identity  for  the  small  z  directed  electric  dipole, 
and  the  reciprocity  relation  (23)  becomes 


R  = 


i  waT. 

Vo 


(44) 


The  power  relations  (30)-(32)  for  the  electrically  small  lossless  z  directed  dipole  simplify 


considerably  to 


ITI2  = 

1  1  Y0  4tt 

Im  ( S )  -  ^  \S\ 


(i  -  |r]2) 

(45) 

,p,2 

T047r  W 

(46) 

V0k2r  * 

Y04*TR  • 

(47) 
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The  reciprocity  relation  (44)  can  be  used  to  eliminate  R  in  (46)  and  (47),  so  that  they  can 
be  written  as 

Im(S)-||S|2  =  ^|T|2  (48) 

T(51  +  lS'HirT''  (49) 

Taking  the  absolute  value  and  then  squaring  both  sides  of  (49)  produces  (48)  after  |T|2  is 
substituted  from  (45)  .  Therefore,  (48)  is  redundant,  and  we  are  left  with  three  indepen¬ 
dent  equations  obtained  from  the  reciprocity  and  power  conservation  relations  for  a  small, 
reciprocal,  lossless  z  directed  electric  dipole  antenna,  namely  equations  (44),  (45),  and  (49). 

Writing  T  in  (44)  and  (49)  as  |T|e'^r  shows  that  the  transmitting,  receiving,  and  scat¬ 
tering  coefficients  can  be  expressed  merely  in  terms  of  the  reflection  coefficient  T  =  |r|e‘^r 
and  the  phase  ipT  of  the  transmitting  coefficient;  specifically 


(50) 

(51) 

S=  ?i  (-|r|ei(2V,r-*r)  +  i). 

(52) 

The  eight  scalars  needed  to  specify  the  magnitudes  and  phases  of  T,  R,  S,  and  T  have  been 
reduced  by  means  of  the  reciprocity  and  power  relations  to  three  scalars  (the  magnitude 
and  phase  of  T  and  the  phase  of  T).  If  the  reference  plane  Ao  in  the  waveguide  feeding  the 
antenna  is  shifted  by  a  distance  £,  the  phases  of  T  and  T  shift  by  /%£  and  2/3o£,  respectively, 
where  o  is  the  propagation  constant  of  the  waveguide  mode.  Therefore,  (52)  confirms  that 
the  phase  of  the  scattering  coefficient  S  does  not  change  by  shifting  the  waveguide  reference 
plane  Ao . 

For  a  lossless  scatterer  rather  than  a  lossless  reciprocal  antenna,  T  =  R  =  0,  |r|  =  1,  and 
(50)-(52)  are  replaced  by  the  one  lossless  power  relation  in  (48)  with  T  =  0: 


|S|2  =  |lm(S) 


or  letting  S  =  |S|e^ 


a  relationship  we  shall  discuss  in  the  next  section. 


(53) 

(54) 


6  INFINITE  LINEAR  PERIODIC  ARRAY  OF  SMALL 
PASSIVE  ELECTRIC  DIPOLE  ANTENNAS  PER¬ 
PENDICULAR  TO  THE  ARRAY  AXIS 

The  transmitting,  receiving,  and  scattering  coefficients  found  in  Section  5  for  a  small,  recip¬ 
rocal,  lossless  z  directed  electric  dipole  antenna  determine  the  propagation  constants  of  the 
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traveling  waves  that  may  exist  on  a  linear  periodic  array  of  these  antennas.  To  obtain  an 
explicit  expression  for  the  dependence  of  the  traveling-wave  propagation  constant  P  on  the 
T,  R,  and  S  of  the  antennas  composing  the  array,  as  well  as  on  the  normalized  spacing  kd 
of  the  antennas,  consider  an  infinite  linear  periodic  array  of  small  z  directed  passive  electric 
dipoles  equispaced  along  the  a:- axis,  each  satisfying  the  scattering-matrix  equations  (41)-(42) 

=  T<  +  REf?  (55) 

bn  =  Ta%  +  SE°n  (56) 

where  the  subscript  n  denotes  the  nth  dipole  in  the  infinite  array.  The  dipole  antennas  are 
both  passive  and  lossless  in  that  they  contain  no  active  or  lossy  material  and  each  of  their 
feed  waveguides  are  terminated  in  a  perfectly  reflecting  load  with  reflection  coefficient  given 
by  Tl  =  e"\  such  that 

aZ  =  (57) 


or 

bZ  =  e-^aZ.  (58) 

Substituting  6ft  from  (58)  into  (55),  solving  for  a£,  then  using  the  resulting  a£  in  (56),  we 
find  bn  in  terms  of  E^1 

bn  =  SeE f  (59) 

where  the  “effective”  scattering  coefficient  is  defined  as 


Se  = 


TR 

e"i -  T 


+  s. 


(60) 


For  a  Unear  array  of  scatterers  rather  than  antennas,  we  have  T  =  R  =  0  and  |r|  =  1,  so 
that  bn  —  SE™  and  Se  =  S.  The  Unear  relationship  between  the  amplitude  bn  of  the  nth 
small  z  directed  dipole  and  the  incident  electric  field  E'p1  at  the  nth  dipole  is  similar  to  the 
linear  relationship  between  the  ampUtude  of  the  electric  dipole  moment  and  the  incident 
electric  field  used  in  the  discrete-dipole  approximation  for  computing  scattering  by  dielectric 
particles  [23],  [24].  In  this  previous  discrete-particle  scattering  work,  a  full  scattering-matrix 
description  was  not  formulated  and  neither  reciprocity  nor  power  conservation  was  explicitly 
applied  to  reduce  the  number  of  variables  and  required  computations. 

The  field  E^1  is  the  electric  field  incident  upon  the  nth  small  dipole  from  all  the  other 
dipoles.  For  an  infinite  array  of  small  z  directed  electric  dipole  antennas  equally  spaced  a 
distance  d  along  the  x  axis,  this  incident  field  E^1  can  therefore  be  expressed  as  a  summation 
of  the  outgoing  fields  of  all  but  the  nth  dipole  antenna  (see  (40)) 


+oo 

*?"  =  £  fc 

j—- OO 


gifcd|j-n 

3  kd\j 


n 


1  + 


+ 


kd\j  —  n|  (kd)2\j  —  n|2J  ’ 


n  =  0,±l,±2,  ...  ±oo.  (61) 


Inserting  E®n  from  (61)  into  (59)  produces  an  infinite  set  of  Unear  equations  for  the  unknown 
outgoing  modal  coefficients 


+oo  piled]  j-n\ 

bn  =  SeYl  bJ  ] 


j—  —  OO 


1  + 


+ 


kd\j-n\  [  kd\j-n\  {kdf\j-n\2\ 


n  =  0,  ±1,  ±2,  ...  ±  oo.  (62) 
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The  effective  scattering  coefficient  Se  defined  in  (60)  has  some  remarkable  properties.  It 
is  the  one  parameter  in  (62),  besides  the  normalized  separation  distance  kd,  that  determines 
the  characteristics  of  the  fields  on  the  linear  infinite  periodic  array  of  small  electric  dipoles. 
Moreover,  if  the  electric  dipoles  composing  the  array  are  reciprocal  and  lossless,  the  T,  R, 
and  S  in  the  equations  (50)-(52)  can  be  inserted  into  (60)  to  show  that  the  magnitude  and 
phase  of  Se  =  |S'e|e1^e  obey  the  simple  relationship 


|*S'e|  =  ^sin^e-  (63) 

Consequently,  ipe  and  kd  are  the  only  two  independent  scalar  parameters  in  the  equations 
(62)  that  determine  the  characteristics  of  the  fields  supported  by  the  infinite  linear  periodic 
array  of  small  z  directed  electric  dipole  antennas.  Because  (63)  demands  that  sini/>e  >  0, 
the  range  of  tpe  is  limited  to 

0  <  V»e  <  7T.  (64) 

We  also  note  that  the  magnitude  and  phase  of  Se  are  not  changed  by  shifting  the  reference 
plane  Ao  in  the  feed  waveguide.  This  result  is  proven  directly  from  (60)  using  the  facts  that 
the  phases  of  T  and  T£  shift  by  2/30t,  the  phases  of  T  and  R  shift  by  fot,  and  the  phase  of 
S  remains  unchanged  (as  explained  in  Section  5)  when  the  reference  plane  Aq  is  shifted  by 
a  distance  £.  (/%  is  the  propagation  constant  of  the  waveguide  mode.) 

We  shall  now  use  (62)  to  determine  an  equation  for  the  traveling-wave  propagation  con¬ 
stant  /?  as  a  function  of  kd  and  ipe.  In  so  doing,  we  shall  also  prove  that  the  relationship  (63) 
is  not  only  a  consequence  of  reciprocity  and  losslessness,  but  it  is  a  relationship  that  must  be 
satisfied  by  the  array  elements  if  the  array  supports  a  traveling  wave.  If  the  array  supports  a 
traveling  wave  with  real  positive  propagation  constant  /?,  the  bj  in  (62)  are  identical  except 
for  a  phase  shift  given  by 

bj  =  boeif3jd  (65) 

which  allows  (62)  to  be  rewritten  in  the  form 


4*oo  ‘\kd\j-n\ 

(kdf  =  S'Y.  ijr^f  L 


j=  —  0O 


(kd)2  +  — 


kd\ 


I  3  ~  n\  |  j  -  n\2\ 


,  n  =  0,  ±1,  ±2,  ...  ±  oo . 


(66) 

Each  equation  in  (66)  for  each  different  n  is  identical  because  the  integer  index  j  in  the 
summation  ranges  from  — oo  to  +00.  Therefore,  simplify  (66)  by  setting  n  =  0  to  get 


(ki)*  =  St-£ 

3=  1 


fei {k+M  ei  (k-/3)dj 


kdi  1 
(kd  2  +  —  - 

3  32 


3  3 

Using  the  summation  formulas  [25,  sec.  1.441,  eqs.  1  and  2;  sec.  1.443,  eqs.  3  and  5] 


(67) 


OO 


E 


cos  na 
n 


5S(1-Lar~ln|2sin(a/2)1’  0<a<2* 


sin  na 


7r  —  a 


0  <  a  < 


(68a) 


n 


2 


(68b) 


^  cos  na  _  7 r2  7r  a  i  a2 

f?x  n2  ~  ~6  ~  ~2  +  7’ 


E 

j=i 


sin  na 


nJ 


7r2a 


7ra2  a3 

T+12’ 


and  the  approximations  (see  Appendix  C) 


0  <  a  <  27r 

0  <  a  <  27T 


(68c) 

(68d) 


V  =  F(a)  «  -0.1381  sin  a  +  0.03212  sin  2 a  -  0.9653a  ln(a/7r),  0  <  a  <  ?r  (69a) 

j=i  n 

F(a)  =  —F(2tt  —  a),  7r  <  a  <  27T  (69b) 

and 

g  cosna  ^  G(a)  _  j  3328  _  0.1424  cos  a +  0.01094  cos  2a  +  0.4902a2  ln(a/7r)  -0.2417a2, 
j=i  n 

0  <  a  <  7r  (69c) 
G(a )  =  G( 27t  —  a),  ir  <  a  <  27t  (69d) 

(67)  is  reduced  to 

(kd)3/Se  =  — ( kd )2  ln[2(cos  fed  —  cos  /?d)]  —  kd  [F(kd  +  (id)  —  F(f3d  —  kd )]  — 

[G(kd  +  (3d)  +  G((3d  -  kd)]  -  i^(fcd)3,  A:d  <  (3d.  (70) 

O 

Writing  Se  as  \Se\ei%i>*,  equating  the  real  and  imaginary  parts  of  (70),  and  dividing  the  two 
resulting  equations,  one  obtains 

\Se\  =  ^sin^e,  kd  <  (3d  (71) 

and 

^(fcd)3  cos  ip e  +  |(— (kd)2  In  [2(cos  kd  —  cos  (3d)]  —  (kd)  [F(kd  +  (3d)  —  F((3d  —  kd)] 

-  [G(kd  +  (3d)  +  G((3d  -  kd)] }  sin  ipe  =  0,  kd<  (3d.  (72) 

We  emphasize  that  the  equation  (71),  which  was  previously  derived  as  (63)  using  the  reci¬ 
procity  and  lossless  power  conservation  relations,  has  been  derived  without  explicit  reference 
to  reciprocity  or  losslessness.  It  emerged  simply,  and  rather  unexpectedly,  as  a  necessary 
condition  for  the  array  equations  (62)  to  have  a  traveling-wave  solution.  For  lossless  scat¬ 
tered  rather  than  lossless  reciprocal  antennas,  Se  =  S  and  (71)  becomes  identical  to  (54), 
which  was  found  previously  from  the  power  relation  for  lossless  scatterers. 

The  equation  (72)  gives  an  implicit  expression  for  the  propagation  constant  (3  of  the 
traveling  wave  supported  by  the  linear  periodic  array  of  lossless,  passive,  electrically  small 
dipoles  as  a  function  of  the  normalized  spacing  kd  of  the  z  directed  short  dipole  antennas 
and  the  phase  ipe  of  the  effective  scattering  coefficient.  Although  a  closed-form  expression 
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cannot  be  found  for  /3d  as  a  function  of  fed  and  tpe  as  was  possible  for  a  linear  periodic  array 
of  lossless,  passive,  acoustically  small  isotropic  antennas  [13],  the  implicit  expression  (72) 
can  be  solved  numerically  for  /3d  given  kd  and  ipe.  Curves  of  /3d  vs  fed  are  plotted  in  Figure  3 
for  different  values  of  tpe  ranging  from  10°  to  170°;  see  (64).  All  the  curves  have  one  branch 
that  begins  at  (3d  equal  to  approximately  1.44,  the  value  for  which  G((3d)  =  0  so  that  (72) 
'  is  satisfied  when  fed  =  0.  All  curves  for  ipe  less  than  about  64°  have  two  branches.  For  the 

curves  for  tpe  less  than  about  58°  one  branch  begins  by  being  extremely  close  to  the  line 
/3d  =  fed.  Of  note  is  the  possibility  of  two  distinct  traveling  waves  that  can  exist  for  the 
same  normalized  spacing  fed  and  effective  scattering  coefficient  phase  tpe.  If.  for  example,  we 
follow  the  line  fed  =  1  we  see  that  it  intersects  the  branches  of  the  ipe  =  45°  curve  in  two 
places  corresponding  to  two  distinct  values  of  the  traveling  wave  propagation  constant  /3. 
The  possibility  that  two  distinct  traveling  waves  can  be  excited  is  quite  different  from  the 
behavior  encountered  when  an  infinite  array  of  isotropic  acoustically  small  transducers  was 
studied  for  which  only  one  traveling  wave  can  be  excited  [13].  By  referring  to  the  definition 
of  Se  in  (60),  we  see  that  the  phase  ipe  (and  magnitude)  of  Se  can  be  adjusted  by  changing 
the  phase  ipL  of  the  reflection  coefficient  of  the  perfectly  reflecting  loads  that  terminate  the 
feed  waveguides  of  the  passive  dipole  elements  of  the  array.  This  adjustment  can  be  done, 
in  principle,  by  changing  the  lengths  of  the  feed  waveguides,  or  by  changing  the  material 
properties  of  the  terminating  loads.  However,  unless  the  dipoles  are  fed  through  a  ground 
plane,  it  may  be  difficult  in  practice  to  appreciably  vary  tpe  in  this  way  and  keep  the  dipole 
spacing  small  enough  to  maintain  a  traveling  wave.  In  general,  the  scattering  phase  ipe  of  an 
electric  dipole  will  change  with  frequency.  In  that  case,  the  propagation  constant  /3  of  the 
traveling  wave  can  be  varied  without  changing  ipe  if  the  separation  distance  d  of  the  array 
elements  is  varied  while  holding  the  frequency  fixed. 

In  a  laboratory,  measurements  of  the  propagation  constant  of  a  traveling  wave  and  its 
variation  with  the  electrical  spacing  of  the  array  elements,  kd,  could  be  efficiently  performed 
by  varying  the  frequency,  keeping  the  dipole  length  and  dipole  separation  constant.  In 
this  case  the  ratio  of  the  length  of  the  dipoles  to  the  separation  of  adjacent  dipoles,  2h/d, 
would  be  fixed  as  would  the  ratio  of  the  wire  radius  to  the  wire  length,  p/h.  Examples 
of  the  theoretical  kd  —  (3d  curves  are  illustrated  in  Figure  4  in  which  three  kd  —  /3d  curves 
corresponding  to  different  values  of  2 h/d  and  with  p/(2h)  =  0.1  are  shown.  The  curves  begin 
at  the  line  /3d  =  fed  and  end  when  (3d  =  n.  A  unique  value  of  the  traveling  wave  propagation 
constant  /3  is  associated  with  a  given  value  of  fed.  The  group  velocity  dfe/d/3  is  positive  for 
all  three  curves.  These  calculations  require  values  of  the  phase,  'tpe,  of  the  effective  scattering 
coefficient  as  a  function  of  h/X  for  fixed  p/h.  These  values  can  be  obtained  using  standard 
scattering  codes  such  as  NEC  [26]  and  CICERO  [27]  since  if  a  thin  short  wire  is  illuminated 
by  a  plane  wave  with  the  electric  field  parallel  to  the  wire,  the  scattered  far  field  is  given  by 
(see  (40))  Eszc  =  belkr/(kr)  with  (see  (59)  )  b  =  SeE ®  so  that  Se  is  the  ratio  of  the  scattered 
field  to  the  incident  field. 

The  theoretical  calculations  of  the  propagation  constant  /3  of  the  traveling  wave  based  on 
(72)  can  be  verified  in  part  by  reference  to  measurements  made  by  Ehrenspeck  and  Poehler 
[1]  who  measured  the  phase  velocity  on  a  Yagi  array  as  a  function  of  the  height  of  the 
directors  for  two  different  spacings  and  two  different  wire  radii.  The  fed  —  /3d  diagram  given 
in  Figure  3  plots  /3d  as  a  function  of  fed  and  tpe,  the  phase  of  the  effective  scattering  coefficient 
Se  =  |  sin  tpe^.  Since  the  phase  velocity  is  equal  to  kd//3d  ,  to  compare  our  theoretical 
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Figure  3:  kd-pd  curves  for  an  infinite  linear  periodic  array  of  short  parallel  electric  dipoles 
perpendicular  to  the  array  axis  with  different  values  of  the  phase  V'e  of  the  effective  scattering 
coefficient. 


20 


M 


Figure  4:  kd-/3d  curves  for  an  infinite  linear  periodic  array  of  short  parallel  electric  dipoles 
perpendicular  to  the  array  axis  with  different  ratios  of  the  dipole  length  to  the  dipole  sepa¬ 
ration  (2 h/d).  The  ratio  of  the  radius  of  the  dipole  to  its  length  is  given  by  p/(2h)  =  0.1. 


21 


results  with  the  experimental  results  of  Ehrenspeck  and  Poehler  we  need  to  obtain  V’e  for 
thin  short  wires  of  different  heights  and  radii.  This  was  done  using  NEC.  In  Figure  5  we  show 
plots  of  ipe  versus  wire  half-length  h/X  for  short  thin  wires  of  three  different  radii.  Figure 
6  shows  curves  of  the  theoretical  phase  velocity  of  the  traveling  wave  on  an  infinite  array 
of  wires  with  spacing  d/X  =  0.2  and  radii  p/X  —  0.012  and  0.024  as  a  function  of  the  wire 
half-length  h/X  along  with  experimental  values  taken  from  Figure  13  of  [1].  Figure  7  shows 
the  corresponding  curves  for  a  spacing  of  d/X  =  0.4.  Considering  the  difficulty  of  making  the 
measurements  the  agreement  is  quite  good,  especially  for  the  shorter  wires.  As  the  length 
of  the  wires  increases  the  assumption  that  the  fields  of  the  wires  are  those  of  elementary 
electric  dipoles  becomes  increasingly  inaccurate. 


7  FINITE  LINEAR  PERIODIC  ARRAY  OF  SMALL 
ELECTRIC  DIPOLES  PERPENDICULAR  TO  THE 
ARRAY  AXIS 


The  traveling-wave  propagation  constant  given  implicitly  by  (72)  was  found  by  solving  ana¬ 
lytically  the  coupled  scattering-matrix  equations  for  an  infinite  linear  periodic  array  of  small 
z  directed  electric  dipoles  equispaced  along  the  x  axis.  We  used  an  infinite  rather  than  a 
finite  number  of  identical  array  elements  to  derive  (72)  because  only  on  an  infinite  array  can 
a  traveling  wave  exist  alone  as  an  exact  homogeneous  solution.  For  a  linear  periodic  array 
with  a  finite  number  of  these  same  antenna  elements,  the  same  traveling  wave  can  exist,  but 
it  is  always  accompanied  by  the  portion  of  the  excitation  fields  that  do  not  couple  to  the 
traveling  wave  and  by  the  fields  generated  through  diffraction  from  both  ends  of  the  finite 
array.  In  this  section  we  consider  such  an  array  of  N  identical  small  z  directed  electric  dipole 
elements  separated  along  the  x  axis  by  a  distance  d,  as  shown  in  Figure  8.  All  the  antenna 
elements  of  this  linear  array  except  the  first  one  are  terminated  in  a  perfectly  reflecting  load 
with  reflection  coefficient  given  by  =  el^L .  The  first  antenna  in  the  array  is  fed  by  an 
incident  waveguide  modal  coefficient  aj. 

We  can  use  the  scattering-matrix  equations  (55)-(56)  to  derive  a  set  of  N  linear  equa¬ 
tions  for  the  N  unknown  outgoing  wave  coefficients  bn  for  each  of  the  array  elements 
(n  =  1,2, 3,  ...  N ).  The  derivation  leading  to  the  equations  in  (62)  for  the  passive  elec¬ 
tric  dipole  elements  in  the  infinite  array  applies  also  to  the  passive  z  directed  electric  dipole 
elements  (n  =  2, 3,  ...  N )  of  the  finite  array.  Thus,  we  can  immediately  write  down  N  —  1 
of  the  required  equations  for  the  finite  array  as 


N  Akd\j-n 
bn-SeJ^h 


1  + 


3kd\j-n\  [  kd\j-n\  (kd)2\j-n\2\ 


=  0,  n  =  2, 3,  . . .  N .  (73) 


The  feed-element  equation  (n  =  1)  can  be  found  by  inserting  E/,1  from  (61)  into  (56)  to  get 


N  Jkd\j-1\ 


1  + 


1) 


kd\j  —  1|  (kd)2\j  -  1|2 


=  Tal 


(74) 


22 


h  IX 


Figure  5:  Phase  ipe  of  the  effective  scattering  coefficient  of  a  thin  perfectly  conducting  wire 
as  a  function  of  the  wire  half-length  hj A. 
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Figure  6:  Theoretical  curves  of  the  relative  phase  velocity  of  a  traveling  wave  on  an  infinite 
linear  periodic  array  of  thin  wires  perpendicular  to  the  array  axis  with  radii  p/ A  =  0.012  and 
0.014  as  a  function  of  the  wire  half-length  h/X  compared  with  experimentally  determined 
values  of  the  relative  phase  velocity  of  a  traveling  wave  on  a  long  Yagi  antenna;  spacing 
d/X  =  0.2. 
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Figure  7:  Theoretical  curves  of  the  relative  phase  velocity  of  a  traveling  wave  on  an  infinite 
linear  periodic  array  of  thin  wires  perpendicular  to  the  array  axis  with  radii  p/X  =  0.012  and 
0.014  as  a  function  of  the  wire  half-length  h/X  compared  with  experimentally  determined 
values  of  the  relative  phase  velocity  of  a  traveling  wave  on  a  long  Yagi  antenna;  spacing  d/ X 
=  0.4. 
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Figure  8:  Finite  linear  periodic  array  of  N  short  electric  dipole  antennas  (represented  by 
the  circles)  perpendicular  to  the  array  axis  separated  by  a  distance  d  along  the  x  axis.  The 
waveguides  of  all  the  elements  of  the  array  are  terminated  in  the  perfectly  reflecting  passive 
load  rL  =  ei^L  except  the  first  element,  which  is  fed  with  an  incident  waveguide  modal 
coefficient  aj. 
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The  solution  to  the  set  of  N  linear  equations  in  (73)-(74)  is  a  trivial  matter  on  a  computer. 
It  is  helpful,  however,  to  specify  the  scattering-matrix  parameters  in  these  equations  in  a 
convenient  form.  First,  we  note  that  it  is  unnecessary  to  specify  the  transmitting  coefficient 
T  on  the  right  hand  side  of  (74)  because  it  is  multiplied  by  the  one  excitation  coefficient  aj, 
which  can  be  chosen  arbitrarily.  Therefore,  for  the  sake  of  simplicity,  let 


(75) 


in  (74).  Secondly,  since  Se  satisfies  (63)  and  (71),  express  Se  in  (73)  as 


Se  = 


0  <  Ipe  <  7T- 


(76) 


Thirdly,  use  (52)  to  rewrite  S  in  (74)  as 


S=!|r|Smael“  +  ?i(l-|r|), 


0  <  a  <  7r,  0  <  T  <  1, 


(77) 


where 


a  =  tan  1 


1  —  cos(2  ij)T  —  ipr) 
sin (2ipT  ~  V’r) 


(78) 


Consequently,  there  are  four  scalar  parameters  ( kd ,  ipe,  a  and  |r|),  in  addition  to  the  number 
of  elements  N  in  the  array,  that  must  be  chosen  to  get  a  numerical  solution  to  (73)-(74). 
We  found,  however,  that  the  numerical  solution  to  (73)-(74)  did  not  change  qualitatively  as 
a  and  |T|  were  varied  over  their  full  range  of  values.  In  other  words,  S  in  (74),  unlike  Se  in 
(73),  is  not  a  critical  parameter  that  must  be  varied  if  one  merely  wants  to  display  different 
representative  numerical  solutions.  In  the  numerical  results  shown  below,  the  values  of  a 
and  |r|  were  set  equal  to  90°  and  .5,  respectively,  so  that  S  has  the  value  of 


S  =  1.125i . 


(79) 


We  also  found  that  the  solution  to  a  multi-element  array  does  not  change  qualitatively  with 
the  number  of  elements.2  The  numerical  results  shown  below  that  do  not  depend  qualitatively 
on  varying  N  are  computed  for  the  representative  value  of  N  —  40.  This  leaves  us  with  just 

two  critical  parameters,  kd  and  ipe,  to  specify  when  solving  for  the  bn  in  equations  (73)-(74). 

Once  the  bn  are  computed  from  (73)-(74),  the  electromagnetic  field  radiated  by  the  array 
can  be  computed  from  the  formulas 

pi kr  N 

E(r)  =r  ~°°  £  6ne-ifcd(n-1)8in9cos^ sin 6  0  (80a) 

KV  _ _ i 

n=l 

H(r)  =r-°°  y0f  x  E(r)  (80b) 

2For  example,  the  number  of  resonant  peaks  in  Figure  9  would  change  in  approximate  proportion  to  the 
number  of  elements  N ,  but  the  height  (~  1.3)  of  the  resonant  peaks  and  their  cut-off  at  kd  ~  2.5  would 
remain  about  the  same. 
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where  9  is  the  spherical  angle  that  the  vector  r  to  the  far  field  makes  with  the  z  axis.  The 
total  power  radiated  by  the  array,  and  the  directivity  function  for  the  array  are  given  in 
terms  of  the  electric  far  field  by 


7 r  Z7T 

£  /  /|rE(r)|SL„stoM*W 


\  47r|rE(r)|2 

D(*)  =  r-Hoo  Tfc? - 

/  /  |ri?(r)|2sin0d0d0 

o  o 


If  only  the  feed  electric  dipole  element  (n  =  1)  were  present  in  free  space,  it  would  radiate 
a  power  equal  to  the  net  power  it  received  from  the  waveguide,  namely 


or  in  view  of  (75)  and  (45) 


fubK i2a  -  in2) 

(83) 

P  _  **Y0 

0  3  k2  ' 

(84) 

,  7 r  27r 

1 1  |rE(r)|r_oosin0d0d0 

(86) 

The  power  ratio 


is  greater  than  or  less  than  1  depending  on  whether  the  array  radiates  more  or  less  power 
than  the  feed  element  radiates  alone  in  free  space.  We  can  also  express  1Z  in  terms  of  the 
magnitude  of  the  in  situ  waveguide-mode  reflection  coefficient  T  i  looking  into  the  feed  electric 
dipole. 

R=idM.  (86) 

where  b\  =  TioJ  and  the  total  power  radiated  can  also  be  expressed  as  P  =  §f?ol°ol2(l  — l^il2), 
the  power  accepted  by  the  feed  element,  because  the  array  is  lossless.  If  1Z  becomes  much  less 
than  1,  the  voltage  required  to  feed  a  fixed  amount  of  power  into  the  array  will  grow  much 
greater  than  that  required  to  feed  the  same  power  into  the  feed  element  when  it  resides  alone. 
In  other  words,  TZ  <C  1  means  that  the  magnitude  of  the  in  situ  waveguide-mode  reflection 
coefficient  |rj|  looking  into  the  feed  electric  dipole  (n  =  1)  is  significantly  closer  to  unity 
than  its  reflection  coefficient  ,|r|  when  all  the  other  dipoles  (n  =  2, 3,  ...  N )  are  removed. 

From  (81),  (82)  and  (85)  together  with  (80a)  the  directivity  and  power  ratio  for  the  finite 
linear  array  along  the  x-axis  of  short  2  directed  dipoles  can  be  calculated  respectively  from 


D(9,4>)  = 


N 

E  Ke 

n=  1 


— ikdn  sin  6  cos  (j)  sin2^ 


TZ  =  SA 
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where 

A  =  ^  53  |6„|2  +  ^  53  53  Re[fen6m]  /  Jo[(n»  -  n)fcdsin  0]  sin3  0d0  (89) 

and  we  have  made  use  of  the  integration  formula  /027r  exp(iu  cos  (f)d(f>  =  2nJ0(u)  with  J0  the 
Bessel  function  of  zero  order.  The  N(N  —  1)  numerical  integrations  required  to  calculate  the 
directivity  and  power  ratio  can  be  avoided  if  instead  of  the  cartesian  xyz  coordinate  system 
we  have  been  using  with  the  x  axis  along  the  array  axis  and  the  z  axis  parallel  to  the  electric 
dipoles  we  switch  to  a  cartesian  x'y'z!  system  with  the  d  axis  along  the  array  axis  and  the  x' 
axis  parallel  to  the  electric  dipoles  (i.e.,  z'  =  x,  x'  =  z,  y'  =  — y).  The  directivity  and  power 
ratio  can  then  be  obtained  in  closed  form: 


and 

where 


D(f,4T) 


^  e— ikdncosO' 

n=l 


(1  —  sin2  0'  cos2  4>') 


2A 


K  =  3A 


(90) 

(91) 


A  =  ^  53  lfcn|2+  53  53  ^Mm*]  (sinc[(m  -  n)kd\  1  - 

^  n=l  n=  1  m=n+ 1  V 


[(m  —  n)kd}2\ 


cos[(m  —  n)fcd]\ 
+  [(m  —  n)kd\ 2  ) 
(92) 


7.1  NUMERICAL  RESULTS  FOR  THE  FINITE  ARRAY 

The  power  ratio  1Z  defined  in  (85)  is  plotted  as  a  function  of  kd  in  Figures  9  and  10  for 
a  40  element  linear  array  of  small  z  directed  electric  dipoles  equispaced  along  the  x-axis. 
In  Figure  9  the  phase  ipe  of  the  effective  scattering  coefficient  has  been  chosen  equal  to  45° 
and  in  Figure  10  it  equals  135°.  In  Figure  9  we  note  that  the  curve  has  a  region  of  large 
oscillations  followed  by  a  region  of  small  moire  equally  spaced  oscillations  and  ending  with 
a  non-oscillatory  portion.  This  seemingly  curious  behavior  can  be  understood  qualitatively 
with  reference  to  the  kd  —  (3d  diagram  of  Figure  3.  For  values  of  kd  between  0.4  and  about 

1.2  two  distinct  traveling  waves  can  be  excited,  one  with  a  phase  constant  close  to  that  of 
free  space,  and  the  other  increasing  from  (3d  equal  to  about  1.44  to  7r.  Apparently  the  two 
traveling  waves  can  both  highly  constructively  and  destructively  interfere  with  each  other 
resulting  in  large  oscillations  of  the  power  ratio.  Also  note  that  one  of  the  traveling  wave 
propagation  constants  is  much  larger  than  k.  The  impedance  mismatch  of  this  traveling 
wave  as  it  encounters  free  space  at  the  ends  of  the  array  will  be  large,  and  an  appreciable 
amount  of  power  will  be  reflected  so  that  a  large  standing  wave  can  be  supported  by  the 
array.  For  kd  between  1.2  and  a  little  more  than  2.6  only  one  traveling  wave  can  be  excited. 
The  propagation  constant  of  this  traveling  wave  is  fairly  close  to  k  so  that  the  impedance 
mismatch  of  this  wave  as  it  encounters  free  space  at  the  array  ends  is  not  very  pronounced 
with  the  result  that  the  corresponding  standing  wave  has  a  relatively  small  amplitude.  For 
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Figure  9:  Power  ratio  versus  kd  for  a  forty  element  linear  periodic  array  of  short  electric 
dipoles  perpendicular  to  the  array  axis  =  45°. 
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Figure  10:  Power  ratio  versus  kd  for  a  forty  element  linear  periodic  array  of  short  electric 
dipoles  perpendicular  to  the  array  axis  with  V’e  =  135°. 
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kd  larger  than  about  2.6  no  traveling  wave  can  be  excited  and  so  there  are  no  oscillations 
of  the  power  ratio.  The  array  radiates  more  or  less  like  the  feed  element  only  in  free  space 
because  none  of  the  feed  element  radiation  excites  a  traveling  wave.  The  scattering  excitation 
of  the  passive  elements  decreases  rapidly  with  increasing  distance  from  the  feed  element. 

In  Figure  10  the  power  ratio  is  plotted  as  a  function  of  kd  for  ipc  =  135°.  Again  the 
kd  —  /3d  diagram  of  Figure  3  can  be  used  to  understand  the  behavior  of  the  curve.  For  kd 
less  than  about  0.7,  although  in  principle  two  traveling  waves  can  exist,  one  has  a  propagation 
phase  constant  almost  exactly  equal  to  that  of  free  space  and  cannot  be  excited  with  a  small  * 

feed  element.  For  kd  between  about  0.7  and  0.8  two  traveling  wave  can  be  excited  resulting 
in  the  more  rapid  oscillations  of  the  power  ratio.  Finally  for  still  larger  values  of  kd  no 
traveling  wave  can  be  excited  and  the  power  ratio  slowly  approaches  the  value  of  1  for  larger 
interelement  spacing. 

An  interesting  feature  of  Figures  9  and  10  is  that  the  resonant  peaks  in  the  power  ratio 
never  get  higher  than  a  value  of  about  1.3.  Therefore,  the  magnitude  of  the  in  situ  waveguide¬ 
mode  reflection  coefficient  |Fi|  given  in  (86)  never  gets  much  less  than  the  magnitude  of  the 
reflection  coefficient  |T|  of  the  feed  element  radiating  alone  in  free  space.  At  kd  values  between 
pronounced  resonant  peaks,  however,  the  power  ratio  becomes  much  less  than  unity,  (Til  gets 
much  closer  to  unity  than  |r|,  and  it  takes  a  much  higher  voltage  to  feed  the  array  with  the 
same  amount  of  accepted  power  as  the  power  accepted  at  the  resonant  peaks  or  when  the 
feed  antenna  is  radiating  alone  in  free  space.  Although  the  curves  in  Figures  9  and  10  were 
computed  for  |r|  =  .5,  the  resonant  peaks  in  the  power  ratio  do  not  exceed  a  value  of  about  2 
as  |r|  — *■  1.  Thus,  operating  at  a  resonant  peak  of  the  array  does  not  eliminate  the  difficulty 
of  feeding  appreciable  power  into  small  electric  dipoles.  To  enable  more  power  to  be  accepted 
and  radiated  by  an  array  of  small  elements,  Veremey  and  Mittra  [10]  have  designed  each 
small  element  to  be  a  resonant  antenna  with  a  reasonably  low  reflection  coefficient.  This 
elemental  resonance  is  superposed  on  the  curves  corresponding  to  those  in  Figures  9  and  10 
for  their  arrays.  When  the  elemental  resonance  occurs  near  an  array  resonance,  the  array 
resonance  appears  to  “split”  into  two  resonances. 

Figure  11  shows  two  far-field  directivity  patterns  computed  from  (82)  for  the  40  element 
array,  one  with  kd  —  2.1  and  i/je  =  30°  and  the  other  for  kd  —  1.9  and  ipe  =  63°.  Since  the 
kd  =  2.1  array  is  a  little  longer  than  the  kd  =1.9  array  it  is  to  be  expected  that  the  directivity 
will  oscillate  more  rapidly  for  the  larger  array  than  for  the  smaller  one.  Also,  referring  to  the 
kd  —  /3d  diagram  of  Figure  3  again,  we  see  that  the  traveling  wave  propagation  constant  for 
the  array  with  kd  =  2.1  and  ipe  =  30°  is  quite  close  to  k  and  so  the  associated  phase  velocity 
is  close  to  the  speed  of  light.  In  contrast,  the  phase  constant  for  the  array  with  kd  =  1.9  and 
%j)e  =  63°  is  considerably  larger  than  k  and  so  the  associated  phase  velocity  is  considerably 
less  than  the  speed  of  light.  The  larger  endfire  directivity  obtained  for  the  array  with  the  • 

larger  phase  velocity  is  an  example  of  a  general  trend  predicted  many  years  ago  by  Hansen 
and  Woodyard  [14]. 

The  difference  in  endfire  directivities  seen  in  Figure  11  can  be  made  much  less  by  changing 
the  number  of  elements  N  for  the  i/)e  =  63°  array  to  a  value  for  which  the  traveling  wave 
launched  from  the  end  of  the  array  and  the  feed  radiation  that  is  not  converted  to  the 
traveling  wave  add  in  phase  (or  nearly  in  phase)  in  the  far  field  at  6  =  0°.  Nonetheless, 
the  maximum  endfire  directivity  attainable  by  varying  the  number  of  elements  in  the  array 
becomes  greater  with  increasing  phase  velocity  of  the  traveling  wave.  We  demonstrate  this 


32 


Figure  11:  Directivity  patterns  for  a  forty  element  linear  periodic  array  of  short  electric 
dipoles  perpendicular  to  the  array  axis  with  kd  =  2.1  and  xbe  —  30°,  and  kd  =  1.9  and 


in  Figure  12  in  which  the  maximum  directivity  in  the  endfire  direction  (6  =  0°),  obtained  by 
varying  N  from  2  to  100,  is  plotted  versus  the  normalized  phase  velocity  of  the  traveling  wave 
[d/c  =  kd/(fid)\  for  three  different  values  of  kd  chosen  with  reference  to  Figure  3  so  that  only 
one  traveling  wave  can  be  excited.  The  general  increase  of  the  maximum  endfire  directivity 
with  increasing  relative  phase  velocity  v/c  is  in  agreement  with  the  prediction  of  Hansen 
and  Woodyard  [14]  (for  closely  spaced  elements,  d/X  <  1/3,  and  ignoring  reflections  of  the 
traveling  wave  from  the  ends  of  the  array)  and  determined  experimentally  by  Ehrenspeck 
and  Poehler  [1].  Of  course,  the  maximum  attainable  endfire  directivity  pf  a  linear  array  of 
N  closely  spaced  isotropic  radiators,  for  which  the  magnitude  and  phase  of  the  excitation 
coefficients  of  the  N  elements  can  be  specified  arbitrarily,  is  much  larger  than  the  maximum 
attainable  endfire  directivity  shown  in  Figure  12  for  the  array  with  one  feed  element  [28]. 


8  LINEAR  PERIODIC  ARRAYS  OF  SMALL  ELEC¬ 
TRIC  DIPOLES  PARALLEL  TO  THE  ARRAY  AXIS 


In  Sections  6  and  7  we  considered  linear  periodic  arrays  of  small  electric  dipole  antennas 
perpendicular  to  the  array  axis.  In  this  section  we  investigate  linear  periodic  arrays  of  small 
electric  dipole  antennas  in  which  the  dipoles  are  aligned  with  the  array  axis.  We  begin  by 
considering  the  possibility  of  traveling  waves  on  an  infinite  array.  We  consider  z  directed 
small  electric  dipoles  and  take  the  z  axis  as  the  array  axis.  The  nth  electric  dipole  in  the 
array  is  excited  by  the  z  component  only  of  the  fields  radiated  by  all  the  other  dipoles  in 
the  array.  Similarly  to  (61),  from  (35a)  and  (39)  the  field  in  the  z  direction  incident  on  the 
nth  dipole  is  given  by 
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Inserting  E°n  from  (93)  into  (59),  assuming  the  existence  of  a  traveling  wave  (65),  and  letting 
n  =  0  yields  the  equation  similar  to  (67) 
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Using  the  summation  formulas  (68c)  and  (68d)  and  the  approximate  summations  (69a)  and 
(69c)  we  obtain  similarly  to  (71)  and  (72), 

|Se|  =  |  sin^e,  kd  <  /3d  (95) 


and 

—~(kd)3  cos  ip e  +  {kd  [ F(kd  +  (3d)  —  F(fid  —  kd)]  + 
o 

[G(kd  +  (3d)  +  G((3d  -  fed)]}  sin  ipe  =  0,  kd  <  fid.  (96) 
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Figure  12:  Maximum  endfire  directivity  versus  relative  phase  velocity  for  a  linear  periodic 
array  of  short  electric  dipoles  perpendicular  to  the  array  axis. 


The  implicit  expression  (96)  can  be  solved  numerically  for  (3d  given  kd  and  ipe.  Curves  of 
(3d,  vs  kd  are  plotted  in  Figure  13  for  different  values  of  ipe  ranging  from  0.1°  to  170°;  see 
(64).  In  contrast  with  the  k/3  -  kd  diagram  of  Figure  3  for  traveling  waves  on  the  infinite 
array  of  short  electric  dipoles  perpendicular  to  the  array  axis,  here  there  is  only  one  branch 
of  the  curve  for  each  value  of  ipe.  All  the  curves  begin  at  (3d  equal  to  approximately  1.44,  the 
value  for  which  G(/?d)  =  0  so  that  (96)  is  satisfied  when  kd  =  0.  Since  a  horizontal  line  for 
a  constant  value  of  kd  intersects  one  of  the  curves  in  only  one  place,  at  most  one  traveling 
wave  can  be  supported  by  an  infinite  array  of  short  electric  dipoles  aligned  with  the  array 
axis. 

Whether  in  fact  it  is  possible  for  a  traveling  wave  to  be  supported  on  a  linear  periodic 
array  of  thin  wire  scatterers  aligned  with  the  array  axis  depends  on  the  actual  value  of  ipe 
for  the  thin  short  wires  of  a  given  length  2 h  and  radius  p  that  compose  the  array.  We  note 
first  that  2 h  <  d  where  d  is  the  separation  of  the  centers  of  adjacent  wires  in  the  array. 
Additionally  the  length  of  the  wires  should  be  considerably  greater  than  the  wire  diameter, 
say  2 h  >  4 p,  for  the  dipole  pattern  (37)  to  hold.  Thus  the  wire  length,  radius,  and  separation 
must  satisfy  the  inequality 

4 kp  <  2 kh  <  kd.  (97) 

Suppose,  for  example,  that  p  =  0.024A.  Then  kp  =  0.15  ,  kh  >  0.30,  and  h/X  >  0.05. 
Referring  to  Figure  5  where  we  have  plotted  the  phase  ipe  of  the  effective  scattering  coefficient 
versus  wire  length  for  wires  of  different  radii  we  can  tabulate  h/ A,  ipe ,  and  2  kh  for  values  of 
h/X  >  0.05 


h/X 

2  kh 

Ipe 

0.05 

0.63 

1.04° 

0.06 

0.75 

1.58° 

0.07 

0.88 

2.26° 

0.08 

1.01 

3.13o 

0.09 

1.13 

4.25° 

0.10 

1.26 

5.66° 

Now  from  the  second  inequality  of  (97)  we  see  by  referring  to  Figure  13  that  the  minimum 
value  of  ipe  necessary  for  a  traveling  wave  to  be  supported  by  an  array  of  wires  of  the  given 
half-length  h/X  is  equal  to  the  value  of  ipe  of  the  k/3  -  kd  curve  that  intersects  the  k/3  =  kd 
line  at  k/3  =  2 kh.  In  all  cases  it  is  found  that  this  minimum  value  of  ipe  is  greater  than  the 
actual  value  of  ipe  tabulated  for  the  wires  of  the  given  length  and  radius.  Since  this  same 
result  is  found  to  hold  for  wires  with  other  values  of  p  we  conclude  that  no  traveling  wave 
can  actually  be  excited  on  an  array  of  short  thin  wires  aligned  with  the  array  axis  (unless 
the  thin  wires  are  loaded  to  substantially  increase  the  value  of  ipe)- 

9  LINEAR  PERIODIC  ARRAYS  OF  SMALL  ELEC¬ 
TRIC  DIPOLES  SKEW  TO  THE  ARRAY  AXIS 

In  this  section  we  consider  linear  periodic  arrays  of  small  electric  dipole  antennas  in  which 
the  dipoles  are  inclined  at  an  angle  60  to  the  array  axis  [29]  (see  Figure  14).  The  electric 
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Figure  13:  kd-(3d  curves  for  an  infinite  linear  periodic  array  of  short  electric  dipoles  aligned 
with  the  array  axis  with  different  values  of  the  phase  ipe  of  the  effective  scattering  coefficient. 
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Figure  14:  Linear  periodic  array  of  short  electric  dipoles  oriented  at  an  angle  6$  with  respect 
to  the  array  axis. 
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dipoles  are  taken  to  be  in  the  z  direction.  To  investigate  the  possibility  of  traveling  waves 
on  an  infinite  array  we  note  that  the  nth  electric  dipole  in  the  array  is  excited  by  the  z 
component  only  of  the  fields  radiated  by  all  the  other  dipoles  in  the  array.  Since  the  z 
component  of  the  electric  field  radiated  by  a  short  z  directed  electric  dipole  is  expressed  in 
terms  of  the  r  and  6  components  of  the  field  by 

Ez  =  Er  cos  6  —  Eg  sin  6,  (98) 


similarly  to  (61),  from  (35a), (35b),  and  (39)  the  field  in  the  z  direction  incident  on  the  nth 
dipole  is  given  by 
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Now  it  is  seen  from  (99)  that  if  sin2  0Q  =  2  cos2  0O  or  60  =  arctan\/2  =ss  54.74°  then  the 
higher-order  distance  terms  vanish  and  we  obtain 
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Inserting  E^1  from  (100)  into  (59),  assuming  the  existence  of  a  traveling  wave  (65),  and 
letting  n  =  0  yields  the  equation  similar  to  (67) 
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Equation  (101)  is  seen  to  be  identical  to  [13,  Eq.(59)],  the  equation  for  a  linear  periodic  array 
of  small  isotropic  electroacoustic  transducers,  with  2/3 Se  instead  of  Se.  Hence  we  obtain 
directly 

|£e|  =  j-smipe,  kd  <  (3d  (102) 


and 


0<i  =  cos'1  (—ie-“"*'fc+ cos  kdj,  kd  <  fid. 


(103) 


Equation  (103)  gives  an  exact  closed-form  expression  for  the  propagation  constant  /3  of 
the  traveling  wave  supported  by  the  linear  periodic  array  of  lossless,  passive,  small  electric 
dipole  antennas  making  an  angle  of  0q  ~  54.74°  with  the  array  axis,  as  a  function  of  the 
normalized  spacing  kd  of  the  dipoles  and  the  phase  ipe  of  the  effective  scattering  coefficient. 
As  kd  —>  0/3d  — >  7r/3  independent  of  ijje.  Figure  15,  taken  from  [13],  contains  curves  of  /3d 
vs.  kd  for  various  values  of  ipe  from  0°  to  179°,  and  shows  that  there  is  at  most  one  traveling 
wave  for  a  given  value  of  and  normalized  spacing  kd.  Referring  to  Figure  5  which  gives 
curves  of  V’e  as  a  function  of  thin  wire  half-length  for  several  different  wire  radii,  we  see  that 
there  are  many  different  short  wires  that  can  be  used  to  form  a  linear  periodic  array  that 
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will  support  a  traveling  wave.  This  conclusion  is  the  opposite  of  that  reached  in  the  previous 
section  when  the  thin  wires  were  aligned  with  the  array  axis.  Although  the  kd  —  /3d  diagram 
for  traveling  waves  supported  by  an  infinite  array  of  short  dipoles  making  an  angle  of  ~  54.74° 
with  the  array  axis  is  identical  to  that  for  traveling  waves  supported  by  an  infinite  array  of 
acoustic  monopoles,  the  power  ratio  and  directivity  for  a  finite  array  of  skewed  electric 
dipoles  are  not  the  same  as  for  an  acoustic  monopole  array  of  the  same  number  of  elements 
because  of  the  polarization  of  the  dipoles.  The  equations  to  obtain  the  finite  array  coefficients 
are  identical  with  [13,  Eqs.(66)-(67)],  that  is  (73)-(74)  with  the  higher-order  distance  terms 
omitted,  and  the  remainder  of  the  analysis  is  the  same  as  for  the  finite  array  of  dipoles 
perpendicular  to  the  array  axis. 


10  CONCLUDING  REMARKS 

A  general  vector  spherical-wave  source  scattering-matrix  description  of  an  electromagnetic 
antenna  has  been  formulated.  Reciprocity  and  power  conservation  are  applied  to  reduce  the 
number  of  independent  variables  in  the  coefficients  of  the  scattering  matrix.  The  general 
formulation  is  then  applied  to  a  short  z  directed  electric  dipole  antenna.  While  the  vector 
spherical-wave  harmonic  source  scattering-matrix  formulation  of  a  general  electromagnetic 
antenna  requires  an  expansion  in  an  infinite  number  of  vector  spherical- wave  harmonics,  the 
electric  field  of  a  short  z  directed  electric  dipole  antenna  can  be  described  with  just  the 
outgoing  Njq  (r)  and  incoming  N^(r)  modes.  The  propagation  constant  of  the  traveling 
wave  on  a  linear  periodic  array  of  short  electric  dipole  antennas  is  shown  to  depend  on  just 
two  parameters  —  the  spacing  (in  wavelengths)  kd  of  the  radiators,  and  the  phase  ipe  of 
the  effective  scattering  coefficient  of  the  radiators.  Three  orientations  of  the  electric  dipoles 
forming  the  array  were  considered:  the  electric  dipoles  perpendicular  to  the  array  axis,  the 
dipoles  aligned  with  the  array  axis,  and  the  dipoles  forming  an  angle  6$  —  arctan  \/2  ~  54.74° 
with  the  array  axis.  A  simple  expression  for  the  propagation  constant  (3  of  the  traveling  wave 
was  found  in  terms  of  kd  and  ipe  in  all  three  cases.  When  the  dipoles  are  perpendicular  to,  or 
aligned  with,  the  array  axis  the  expression  for  the  propagation  constant  is  a  transcendental 
equation  which  can  be  easily  solved  numerically.  When  the  dipoles  are  oriented  at  an  angle 
of  do  =  arctan  y/2  with  the  array  axis,  a  closed  form  expression  for  the  propagation  constant 
can  be  obtained  that  is  identical  to  the  expression  that  holds  for  the  traveling  wave  on  an 
array  of  acoustic  monopoles.  Curves  of  the  functional  dependence  of  /3d  on  kd  for  different 
values  of  the  phase  of  the  effective  scattering  coefficient  i/)e  are  given  for  all  three  cases. 
For  the  perpendicular  array,  the  relative  phase  velocity  of  the  traveling  wave  calculated 
from  the  kd  —  /3d  diagram  using  numerical  values  for  the  phase  of  the  scattering  coefficient 
for  thin  short  wires  obtained  with  the  NEC  scattering  code  are  in  close  agreement  with 
measurements  made  by  Ehrenspeck  and  Poehler  [1].  If  the  electric  dipoles  aligned  with 
the  array  axis  are  simple  short  thin  wires,  the  phases  of  the  effective  scattering  coefficients 
obtained  from  the  NEC  code  taken  in  combination  with  the  corresponding  kd  —  (3d  diagram 
preclude  the  possibility  any  traveling  waves  being  excited.  For  a  Yagi-Uda  finite  array  of  N 
electric  dipoles  perpendicular  to  the  array  axis  the  two  parameters  kd  and  ipe  are  the  only 
critical  variables  in  the  solution  for  the  radiation  coefficients. 
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Although  the  scattering-matrix  analysis  of  arrays  developed  in  this  paper  was  applied 
only  to  arrays  of  short  electric  dipoles,  the  basic  formulation  applies  to  any  array  composed 
of  linear,  reciprocal,  lossless  array  elements.  A  preliminary  investigation  suggests  that  an 
extension  to  arrays  consisting  of  small  spheres  of  a  given  permittivity  and  permeability  is 
feasible.  Such  spheres  can  be  described  in  terms  of  a  pair  of  crossed  electric  and  magnetic 
dipoles  and  would  require  four  vector  spherical  wave  function  modes  (two  outgoing  and  two 
incoming)  for  their  description.  An  extension  to  electrically  larger  array  elements  would 
require  more  spherical  modes  to  represent  their  radiated  and  scattered  fields,  and  thus  more 
sophisticated  procedures  to  obtain  an  efficient  numerical  solution  for  both  infinite  and  finite 
arrays.  In  addition,  it  seems  likely  that  the  present  analysis  could  be  extended  straight¬ 
forwardly  to  two-dimensional  linear  arrays  [30]  and  to  arrays  that  do  not  lie  in  a  straight 
line  such  as  closed-loop  traveling- wave  arrays  [12],  and  to  surface  or  volumetric  arrays  of 
acoustic  and  electromagnetic  radiators  and  scatterers  (including  small  resonant  scatterers) 
[10],  [31]— [33]. 


A  TRAVELING  WAVE  INEQUALITIES 

Even  though  the  inequalities  in  (1)  are  stated  in  a  number  of  papers,  texts,  and  antenna 
handbooks,  we  have  not  seen  a  proof  of  (1)  for  a  general  periodic  array  of  electrically  small 
elements.  The  second  inequality  in  (1),  namely,  (3d  <  7r,  follows  from  the  fact  that  as  /3d 
becomes  greater  than  7 r  the  traveling  wave  can  merely  be  re-expressed  as  another  traveling 
wave  propagating  in  the  negative  z  direction  with  a  f3'd  <  7r.  To  see  this,  simply  note  that  a 
periodic  function  that  satisfies  p(x,  y,z  +  d)  =  p(x,  y,  z)el^d  with  7r  <  [3d  <  2tt  also  satisfies 
p(x,  y,z  +  d )  =  p(x,  y,  z)e-',f3'd  with  0  <  (3'd  =  2ir  -  [3d  <  7r,  which  is  a  wave  traveling  in 
the  — z  direction.  For  2i r  <  fid  <  Sir,  the  traveling  wave  can  again  be  written  as  a  wave 
propagating  in  the  +z  direction  with  a  propagation  constant  having  a  value  less  than  7 r.  For 
<  f3d  <  4ir,  the  traveling  wave  reverts  to  a  wave  propagating  in  the  —z  direction  with  a 
propagation  constant  having  a  value  less  than  7 r,  and  so  on.  Therefore,  the  full  range  of  the 
absolute  value  of  the  propagation  constant  of  any  traveling  wave  can  always  be  confined  to 
a  value  equal  to  or  less  than  ir  (see  Footnote  1). 

The  first  inequality  in  (1)  implies  that  the  traveling  wave  is  a  slow  wave  in  that  the  phase 
velocity  of  the  traveling  wave  is  equal  to  or  less  than  the  speed  of  fight.  This  first  inequality 
can  be  proven  for  a  general  linear  periodic  array  of  lossless  passive  electrically  small  elements 
by  looking  at  the  radiation  pattern  of  such  an  array  as  the  number  of  elements  in  the  array 
becomes  infinite.  The  far-zone  electric  field  E(r),  radiated  by  a  linear  array  consisting  of 
an  infinite  number  of  identical  antenna  elements  separated  by  a  distance  d  from  z  =  —00 
to  z  =  +00,  such  that  each  successive  radiating  element  has  the  same  complex  far-zone 
electric  field  pattern  2^(0,  <j>)  (which  is  proportional  to  the  “embedded”  element  pattern) 
phase  shifted  by  an  amount  e,/W,  is  given  by 
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measured  from  the  positive  z  axis.  The  far-field  formula 
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used  in  (104)  can  be  derived  under  the  assumption  that  the  sources  are  either  confined  to 
a  finite  region  or  decay  exponentially  at  large  distances  from  the  origin.  Thus,  the  small 
positive  imaginary  parts,  i q1  and  iq2,  of  (3  have  been  inserted  into  (104)  to  give  the  sources 
an  exponential  decay  as  z  approaches  plus  and  minus  infinity,  respectively.  At  the  end  of 
the  derivation,  q1  and  q2  are  allowed  to  approach  zero. 

Use  of  the  summation  formula,  Y,m=o  xTn  =  1/(1  —  x),  |ar|  <  1,  converts  (104)  to 
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In  the  limit  as  q1  and  q2  approach  zero,  the  quantity  inside  the  square  brackets  of  (105)  equals 
0  for  (/?  —  kcos0)d  ^  2n7r,  and  equals  oo  for  ((3  —  k  cos  6)d  =  2mr,  n  =  0,  ±1,  ±2,  ...  ±  oo. 
However,  the  integration  of  the  bracketed  quantity  in  (105)  with  respect  to  the  variable 
u  =  (/?  —  k  cos  9)d,  over  limits  that  cover  just  the  singularity  at  u  —  0,  is  equal  to  47r  as  qx 
and  q2  approach  zero.  The  same  is  true  for  limits  of  u  that  just  cover  the  other  singularities 
at  2mr,  n  =  ±1,  ±2,  ...  ±oo.  Therefore,  as  ql  and  q2  approach  zero,  the  bracketed  quantity 
in  (105)  is  simply  an  infinite  sum  of  delta  functions  given  by 


Akr  oo 

E(r)  r  ~°°  47r — F(d,  4>)  5[((3  —  k  cos  6)d —  2mr].  (106) 

r  n=— oo 


Each  of  the  delta  functions  in  (106)  represent  power  radiated  (actually  an  infinite  amount 
of  power  radiated  because  the  integral  of  the  square  of  a  delta  function  is  infinite)  by  the 
array  at  the  angles  9  given  by  the  zeroes  of  the  arguments  of  the  delta  functions,  provided 
the  far-field  pattern  (F{9,  <f>)  does  not  have  a  null  at  these  9  angles  for  all  values  of  0.  The  far- 
field  pattern  of  small  (nonresonant)  antenna  elements  will  contain  no  such  nulls.  Therefore, 
since  a  passive  array  carrying  a  lossless  traveling  wave  cannot  radiate  and  still  maintain  the 
traveling  wave,  for  arrays  of  electrically  small  antenna  elements,  the  arguments  of  the  delta 
functions  must  not  be  zero  for  any  9  other  than  9  =  0,  the  direction  of  propagation  of  the 
traveling  wave.  In  particular,  if  the  argument  of  the  n  =  0  delta  function  cannot  be  zero  for 
all  9^0,  then  (3  >  k]  that  is,  the  traveling  wave  is  a  “slow”  wave  because  its  phase  velocity 
is  less  than  or  equal  to  the  speed  of  fight.  Also,  (106)  shows  that  there  can  be  no  traveling 
wave  unless  kd  <  n  even  if  /3d  is  allowed  to  range  from  — oo  to  +oo,  provided  again  that 
there  is  no  null  in  T{9,  0)  at  9  =  9o  >  0  for  all  <j>.  If  there  is  such  a  null,  then  a  fast  wave 
with  f3d  =  kd  cos  9o  may  exist  even  if  kd  >  ir. 

In  all  we  have  proven  that 

kd  <  (3d  <  7r  (107) 

for  a  general  infinite  linear  periodic  array  of  lossless  passive  antenna  elements  for  which 
there  is  no  0  >  0  such  that  T{9, 4>)  =  0  for  all  <j>  —  in  particular,  for  electrically  small 
(nonresonant)  antenna  elements.  An  analogous  proof  shows  that  the  same  inequalities  apply 
also  to  the  propagation  constant  of  the  traveling  waves  on  linear  periodic  arrays  of  lossless 
passive  transducer  elements  with  no  9  >  0  such  that  the  far-field  element  pattern  T{9, 4>)  —  0 
for  all  </>. 

If  the  element  pattern  T(9,  </>)  of  each  antenna  element  does  have  a  null  at  some  9  —  9o  >  0 
for  all  0,  presumably  a  fast  traveling  wave  with  phase  velocity  greater  than  the  speed  of  fight 
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(/?  =  k  cos  $o  <  k )  that  does  not  radiate  can  propagate  along  the  array.  If  such  an  array 
were  used  as  a  scanning  phased  array  antenna,  the  array  antenna  would  have  a  “blind  spot 
at  the  scan  angle  equal  to  the  null  angle  6q  of  the  element  pattern  [34]— [37]. 

As  a  conclusion  to  this  appendix,  we  shall  give  an  alternative  proof  that  j3  >  k  for  far- 
field  element  patterns  having  no  nulls  with  respect  to  9  >  0.  This  alternative  proof,  which 
also  reveals  that  traveling  waves  on  a  linear  infinite  periodic  array  must  decay  exponentially 
with  radial  distance  from  the  axis  of  the  array,  begins  by  expanding  the  fields  outside  the 
elements  of  the  array  in  terms  of  cylindrical  waves.  To  avoid  using  vector  cylindrical  waves, 
assume  an  array  of  acoustic  elements  and  express  the  acoustic  pressure  outside  the  elements 
carrying  a  traveling  wave  in  a  complete  set  of  outgoing  scalar  cylindrical  waves  [22,  ch.  6] 

+O0  +°° 

p(p,<M=  E  /  CnWPiKp^e'Td'y  (108) 

n=— oo__„ 


where  the  H^(K,p)  are  the  cylindrical  Hankel  functions  of  the  first  kind  and 

K,  =  (k2  —  72)* 


(109) 


is  chosen  positive  real  or  positive  imaginary. 

Applying  the  periodic  boundary  condition  for  the  traveling  wave 

p(p,</),z  +  d)=p(p,<j),z)e'pd  (HO) 

and  using  orthogonality  of  the  em4>  and  e<7*  functions  converts  (108)  to 


+oo  +0O 


(111) 


ri——oo  m=— oo 


with 


Km  =  [k2  -  ((3 +  2wm/d)2]2.  (112) 

The  cylindrical  modes  in  (111)  are  also  the  “Floquet  modes”  in  the  parallel  plate  waveguides 
of  width  d  satisfying  the  periodic  boundary  condition  (110). 

The  cylindrical  far-field  pattern  (that  is,  the  far-field  pattern  obtained  by  letting  the 
length  of  the  array  approach  infinity  in  the  ±z  directions  before  letting  p  — >  oo)  can  be 
found  by  replacing  the  Hankel  functions  in  (111)  by  their  large  argument  asymptotic  forms 
to  get 

+oo 

1  (113) 


in  which 


I  9  +0° 

V  7V^m  n=—oo 


The  cylindrical  far-field  pattern  in  (113)  can  be  rewritten  simply  as 


4-00 


(114) 


(115) 


J71——00 


<* 


J/ 
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where 


(116) 


t 


Offi  —  tan 


Km 

(3  +  27T  m/d 


If  9m  is  real,  the  mth  mode  in  (115)  represents  a  cylindrical  wave  propagating  along 
a  cone  with  half  angle  equal  to  9m.  Thus,  for  (115)  to  represent  the  cylindrical  far-held 
pattern  of  a  nonradiating  traveling  wave,  all  the  9m  must  be  imaginary.  That  is,  the  K,m 
must  be  imaginary  for  all  nonzero  Cm(<f>).  Now  (115)  also  reveals  that  Cm((f> )  cannot  be 
zero  for  all  <f>  unless  p(r,  0,  <j>)  (defined  by  its  analytic  continuation  if  9  is  complex)  has  no 
discrete  cylindrical  wave  traveling  in  the  direction  0  =  6m,  or,  equivalently,  unless  the  far- 
held  pattern  of  each  element  in  the  infinite  array  has  a  null  at  6  =  6m  (for  all  <f>).  In  any  case, 
the  cylindrical  far  held  of  a  nonradiating  traveling  wave  with  el/3z  dependence  on  a  linear 
infinite  periodic  array  decays  exponentially  with  radial  distance  as 


exp  i (3z  —  yfi/3~+  2-Kmi/d)2  —  k2  p 

VP 


(117) 


as  p  — ►  oo,  where  the  integer  mi  is  the  index  of  the  slowest  exponentially  decaying  Floquet 
mode  that  is  excited  by  the  array  elements. 

Suppose  the  element  far-held  pattern  has  no  nulls  for  real  6  (and  all  (j>)  and  Co  =  0.  Then 
if  (3  <  k,  there  must  be  a  null  in  the  far-held  pattern  for  real  9  (and  all  <f>)  since  Co  =  0. 
Consequently,  if  C0  =  0  and  the  element  far-held  pattern  has  no  nulls  for  real  9  (and  all  <fr), 
then  (3  >  k  and  the  only  way  that  Co  can  be  zero  is  for  there  to  be  a  null  at  a  complex  value 
of  9  (and  all  <f>)  in  the  analytically  continued  element  far-held  pattern. 


B  VECTOR  SPHERICAL  WAVE  FUNCTIONS 

In  this  Appendix  we  give  the  definitions  of  the  vector  spherical  wave  functions  M  and  N  used 
in  this  report.  The  definitions  are  those  used  by  Billy  Brock  [16]  adapted  for  the  exp(—iu>t) 
time  dependence  used  here: 


M£(r)  =  U 


(21  +  !)(/  -  m)! 
^  4ttI(1  +  1  )(l  +  to)! 


-—  zY i)  (hr)P{n(  cos  0)eim*0 
sin  0 


+  smfe»(*r)-/rM 


z—cos  6 


3im*0 


n£m  =  i. 


(21  +  1)(Z  —  to)! 

^  4wl(l  +  1)(Z  +  m)\ 


-j— *(*  +  l)F>im(cos  0)eim^? 

KT 


j__a 

kr  dr 

J_d_ 
kr  dr 


|^rz^(A;r)j  sin9-^-P[n(x) 


^m<f>Q 


x=cos  0 


+Tr7Tr  KH  C08«)e«^, 


(118) 


>  • 


(119) 
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In  (118)  and  (119)  2,(1)  and  z,(2)  are  the  spherical  Bessel  and  Hankel  functions  ji  and  fy(1), 
respectively,  and  Pzm(x)  is  the  associated  Legendre  function  given  by 


ri" 

/?•  <x)  =  (-ir(i  >  o, 

P,-W(x)  =  >  °. 


where  p(x)  is  the  Legendre  function 


(120a) 


(120b) 


(121) 


The  vector  spherical  wave  functions  M  and  N  given  by  (118)  and  (119)  can  be  defined 
in  terms  of  the  normalized  radially-independent  vector  spherical  harmonic  function  X;m  of 
Jackson  [15]  by 

z^\kr)Xlm(0,<l>)  (122) 

and  1 

Ng(r)  =  ±V  x  *«(*r)Xfcn(M),  (123) 

and  are  related  to  the  and  vector  spherical  wave  functions  of  Hansen  [20]  by 

M«(r)  =  iFii(r)  (124) 

and 

Ng(r)  =  >F&(r)(f)  (125) 

where  the  superscript  i  =  1  when  c  =  1  and  the  superscript  i  =  2  when  c  =  3.3 


C  SUMMATIONS  OF  TRIGONOMETRIC  SERIES 

In  this  Appendix  we  discuss  the  approximations  used  for  the  sums  of  the  trigonometric  series 


^  sin na  _ 

F(°)  =  E-^~>0<a<7r> 

j=i  n 


(126a) 


.  .  v COS  Tld 

G(°)  =  E-3->0<a<7r- 


(126b) 


Closed  form  expressions  are  not  available  for  these  sums.  The  IMSL  least-squares  approxi¬ 
mation  program  FNLSQ  was  used  to  compute  the  approximations 


F(a)  -0.1381  sin  a  +  0.03212  sin  2a  -  0.9653a  ln(a/7r),  0  <  a  <  tt, 


(127a) 


3These  relations  between  the  vector  spherical  harmonics  of  Brock  and  those  of  Hansen  are  given  incorrectly 
in  the  original  form  of  [16]  but  have  been  corrected  by  Brock  in  an  errata. 
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and 


G(a)  «  1.3328  —  0.1424 cosa  +  0.01094 cos 2a  +  0.4902a2 ln(a/7r)  —  0.2417a2,  0  <  a  <  ir. 

(127b) 

Figures  16  and  17  show  F(a)  and  G(a),  respectively,  calculated  with  1000  terms,  along  with 
their  least  squares  approximations.  The  agreement  of  the  approximate  with  the  exact  curves 
is  excellent.  It  will  be  noted  that  although  G'(a)  =  —F(a),  when  the  approximation  for 
G(a)  is  differentiated  the  result  differs  slightly  from  the  negative  of  the  derivative  of  the 
approximation  for  F(a).  While  it  is  possible  to  obtain  an  approximation  for  F(a)  by  taking 
the  negative  of  the  derivative  of  the  approximation  for  G(a),  the  direct  least-squares  fit 
approximation  for  F(a)  that  we  have  used  distributes  the  errors  in  the  approximation  more 
uniformly  over  the  interval  from  0  to  i r. 
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